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ABSTRACT 

Models of poststarburst (or "K+A" ) galaxies are constructed by combining fully three-dimensional 
hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. 
Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated 
as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, 
orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that 
the lifetime and properties of the K-f A phase are strong functions of the merger scenario. K-l-A 
durations are generally < 0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which 
is obtained only in rare cases, owing to a wide variation in star formation histories resulting from 
different orbital and progenitor configurations. Combined with empirical merger rates, the model 
lifetimes predict rapidly-rising K-l-A fractions as a function of redshift that are consistent with results 
of large spectroscopic surveys, resolving tension between the observed K-l-A abundance and that 
predicted when one assumes the K+A duration is the lifetime of A stars (^1 Gyr). These simulated 
spectra are spatially resolved on scales of about 1 kpc, and indicate that a centrally-concentrated 
starburst causes the Balmer absorption strengths to increase towards the central few kiloparsecs 
of the remnant. The effects of dust attenuation, viewing angle, and aperture bias on our models are 
analyzed. In some cases, the K+A features are longer-lived and more pronounced when AGN feedback 
removes dust from the center, uncovering the young stars formed during the burst. In this picture, 
the K+A phase begins during or shortly after the bright starburst/ AGN phase in violent mergers, 
and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the 
gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates 
of K+A incidence as a function of merger scenario. 

Subject headings: galaxies:evolution, galaxies interactions, galaxies:starburst, methods: numerical 



1. INTRODUCTION 

Generally speaking, bright galaxies fall into two broad 
visible-spectrum categories: star-forming, blue spiral 
galaxies, and passively-evolving, red elliptical galaxies 
(|Hubblelll926[ ). Two long-standing questions for study- 
ing galaxy formation are: How does each type of galaxy 
form? Do galaxies transform from one type to another, 
and if so, how? 

Early studies of the relationship between these classes 
as a furi c tion o f loo kback time focused on g alaxy clusters. 
lOemleii ()1974[ ) and lButcher fc Oemleii (|l97S) concluded 
that intermediate-redshift (z ~ 0.5) clusters have a much 
larger fraction of blue, star-forming galaxies t han nearby 
cluste rs: the "Butcher-Oemler effect" . Dress ier fc GunnI 
(jl9.83!) found that in one cluster, these blue cluster mem- 
bers were not "normal" spiral galaxies as we might see in 
the local Universe, but had spectra indicating either nu- 
clear activity (e.g. Seyfert) or a recently quenched star- 
forming episode. The latter, so-called "E+A" galaxies, 
named as the combination of an early-type galaxy spec- 
trum (E) with an A type stellar spectrum, were fitted 
best by an old stellar population plus a starburst com- 
ponent whose development ended within the last ^ 1 
Gyr; they have strong Balmer absorption features but 
weak or undetected nebular emission lines. Also known 
as "K+A" (for K type stars) or "poststarburst" galaxies, 
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E+A galaxies were seen as occupying an important tran- 
sitional state between the star-forming and passive popu- 
lations. These objects are commonly selected from spec- 
troscopic surveys as having H(5 equivalent width greater 
than some value (typically roughly SA) in absorption and 
[O II] emission equivalent width less than roughly 2A. 
The latter limit corresponds to a cut in specific star for- 
mation rate at roughly 10~"'^^yr~^, or a star formation 
rate below about 2 Moyr"^ for the most ma ssive galax- 
ies we consider in this work (jKennicu t^ [T99l . 

Numerous studies of K+A galaxies have been un- 
dertaken to explore this paradigm of galaxy evolu- 
tion. K+A galaxies were found to be a significant 
fraction (^ 0.2) of rich cluster members at 0.2 < 
z < 0.6 (Couch & Sharpies 1987; 'Dressier & Gunn"199Q|: 
f^brican t ct al. 19 91; BcUoni et al. 1995; Balogh ct al] 
19991: iPoggianti et al.lll999D . but are fewer t han 1% of 
cluster members at z < 0.1 (|DressleH 119871 ). Thus, a 
cluster-specific mechanism was thought to be responsi- 
ble for quenching the star formation in these galaxies at 
intermediate redshifts, leaving fewer potential progeni- 
tors for the present time. Unbiased spectroscopic sur- 
veys revealed that although K+As seem to be a sizable 
fraction of cluster galaxies, they also e xist in substan- 
tial n umbers in l ower density regions (iZabludqff et ahl 
1991 [B lake et al|[200l iGotol [2007ai lYan et al.l l2009aF 
Vergani et al.i 12010( 1 out to z ~ 1. In fact, the ma- 
jority of K+A g alaxies in the local u niverse do not re- 
side in clusters (|Quintero et al.l [20041) . Thus, a cluster- 
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specific mechanism ca nnot be the only way to make a 
K+A galaxy (also see lYan et al.ll20d9al ). and a system- 
atic study of alternative mechanisms is warranted. 

Detailed studies of the morphology, kinematics, and 
metal abundances of nearby poststarburst galaxies (pri- 
marily in the field) indicate that many of t hem are likely 
proge n itors of bul g e-doin ii iated galaxies (iNortqn et al. 
2OOII: iBlake et al.| [200I iQuintero et al.l 12004 iGoto 



2007bHYang et al.ll2008ir Thus, studying K+A galax- 



ies may give clues to the evolution of galaxies from 
star-forming disks into quiescent ellipticals in a va- 
riety of enviroii i nents at different cosmological times. 
iLaverv fc Henrvl (|1988[ ) proposed galaxy-galaxy interac- 
tions as a mechanism for driving K-f-A formation. A 
number of previously-mentioned studies of nearby K-|-As 
point out that a substantial fraction exhib it tidal features 
iZabludoff et al."1996; Chang et al." 2001) or companion 
galaxies (Yama uchi e t al. 2008), suggesting that galaxy 
interactions or mergers indeed play a role in K-l-A forma- 
tion. In addition, hydrodyn amic simulations (with sim- 
ple s tar-formation laws, e.g. iMihos fc HernquistI [l994bl 
Il996( ) of gas-rich galaxy mergers can lead to global star 
formation histories (SFH) characterized by strong, short 
bursts. Using simple stellar population synthesis (SPS), 
these bursts can lead to spectra with the lack of nebular 
star-formation indicators and strong Balmer absorption 
and that are found in poststarburst galaxies. 

These studies and others led to proposed K-l-A forma- 
tion mechanisms that fall broadly into two classes: (1) 
interaction of a star-forming disk with another galaxy, 
and (2) processes specific to massive clusters and ha- 
los. While a non-cluster-specific mechanism is required 
for K+As in the field, the relative inefficiency of ma- 
jor mergers in systems of galaxies with high relative ve- 
locities implies that galaxy-galaxy interactions may not 
currently be sufficient for creating the large number of 
K-t-As in the highest-mass halos. 

Ram-pressure stripping ()Gunn fc GottlUlTl . the re- 
moval of gas from disks traveling at high-speeds through 
a hot intracluster medium (ICM), is a likely candidate 
for reducing or ending the star formation in galaxies 
that fall into clusters. Although this process alone may 
not be violent enough to trigger a K-|-A phase, numer- 
ical work indicates that such galaxy-cluster interactions 
may lead to compression of the disk gas, offsets of the 
disk gas with respect to the galaxy's central potential, 
or infall of st ripped gas that remained bound to the sub- 
halo (e.g. seelSchulz fc Struckl [200TI: IVollmer et ahllWiTI : 
Ivan GorkomI \2mW ~ These may lead to enhancements 
of the star formation rate, and thus a better chance 
of seeing s trong Balmer lines in a galaxy's spectrum 
(jPressler fc Gunn 1983). However, these systems can- 
not dominate the population of bright K+A galaxies, 
because the vast majority o f galaxies at thes e masses 
(> IO^Mq) are centrals fe.g. lYang eralll2009[ ). so they 
undergo their transition from star-forming to passive as 
central galaxies. 

In addition, for galaxies in halos above ^ 10^'^ Mq, 
cosmological accretion may lead to virial shocks 
and heating (jRees fc Ostrikeri Il977| ) that shut down 
star formatio n. Numerical simulations of this phe- 
nomenon (e.g. IKeres et al.ll2005l:IDekel fc Birnboiml[2006l: 
iBirnboim et al.ll2007D indicate that galaxies in massive 
halos may experience a powerful burst-quench cycle lead- 



ing to a post-starburst phase. However, this mode is 
likely not a dominant contributor to the K-|-A popula- 
tion for several reasons. In particula r, gas will continue 
to co ol and rain down on the galaxy (jKeres fc Hernquis^ 
'2009") , so it is not certain that a complete (or fast enough) 
cessation of star formation will occur. In addition, it is 
not clear if this mode can drastically alter galaxy mor- 
phology in a manner consistent with field K-l-A galaxies. 

In this work, we focus on K-l-A formation via galaxy- 
galaxy interactions. Mergers have long been thought to 
be responsible for destroying stellar disks and turning 
them into di spersion-domin ated ellipticals: the "merger 
hypothesis" (jToomrd I1977D . However, recent advances 
in numerical simulation methods have allowed the study 
of galaxy evolution via mergers in a significantly more 
quantitative way. For example, it has recently become 
possible to "age-d ate" the t ime since the starburst in 
K-l-A galaxies (e.g. lYang et al. 2008). and merger simu- 
lations can be used to test these methods by including a 
more sophisticated treatment of dust attenuation. 

In particular, using physical models also allows the 
study of realistic SFH outcomes for a given set of merg- 
ers. Empirical models used in studies of K-f A galaxies 
have often assumed very simple forms for the SFH in- 
duced by a merger. These SFHs imply "K-l-A lifetimes" 
of order the lifetime of an A type star, 0.5-1.5 Gyr. While 
it is true that for a given galaxy, the merger-induced 
burst can mimic this trivial single-aged stellar popula- 
tion (SSP), it is not clear if this is true for the "typical" 
merger involving a given disk. In fact, with fixed pro- 
genitors, variations in the orbital parameters of a merger 
can lead t o an enormous variety of dynam ical and SFH 
outcomes (|Barneslll99aiBarnes fc Herngu ist 1996). The 
merger or event rate inferred by surveys of K+A galax- 
ies are sensitive to this distribution of the K+A duration 
timescale. In the present work, we quantify this variation 
on the inferred K+A lifetimes of simulated mergers. 

In addition, the mechanism for the rapid shutoff of 
star formation is not known precisely. Feedback from 
star formation and active galactic nuclei (AGN) activity, 
presumably leading to the expulsion of the remaining gas 
reservoir, has been inferred indirectly in a wide var i ety of 
observations (iKavirai et al.l 120071: iTremonti et al.l 120071: 
iBrownet ani2009D . By varying the feedback prescrip- 
tions and strengths in the numerical models, simulations 
can test a large variety of these processes in a systematic 
way. More generally, by studying poststarburst galaxies 
in the context of the extragalactic zoo and standard cos- 
mology, we hope to study their implications for the role 
of mergers in driving galaxy evolution. 

In this paper, we define and study a set of K+A galaxy 
models which are the remnants of simulated galaxy merg- 
ers. We have seen that the characteristic rapid inflow 
and consumption of gas during a major gas-rich merger 
is likely to be the cause o f some K+A galaxies, and a 
recent studv bvl^yild et all (|2009[) used simulations from 
[Johansson et alF i 20091 ) to perform a similar comparison 
of the models to poststarburst data; the reader is en- 
couraged to see those papers for another detailed inves- 
tigation of the characteristics of poststarburst galaxies in 
merger simulations. Our co nclusions are broa dly consis- 
tent with those reported in iWild et all (|2009D . 

The present work applies different spectrum synthe- 
sis and selection methods to a broad set of simulations 
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bv lCox et al.l ()2006al ) to systematically quantify the vari- 
ation of K+A properties with different merger parame- 
ters. In particular, we utilize three-dimensional radiative 
transfer to study the post-starburst stellar populations 
using standard line strength indicators. This procedure 
allows us to arbitrarily vary observing parameters such 
as the viewing angle and spectrograph fiber size and re- 
alistically measure the line strengths for a wide range of 
simulated mergers that include processes such as star for- 
mation, feedback, accretion onto a central supermassive 
black hole, and dust attenuation. 

The primary goal of this work is to reconcile observed 
K-l-A numbers with realistic merger models and well- 
motivated merger rates. In fj2l we outline in detail our 
methodology for generating spectral line catalogs from 
numerical simulations. Readers wishing to bypass these 
details can skip to the relevant parts of S|3l 21 ^3 oi' iHl 
In Sj3l we study the basic ramifications of our methods, 
highlighting the resulting variations owing to viewing an- 
gle, aperture bias, the presence of a bulge, and AGN feed- 
back. In ^ we systematically calculate the lifetime of 
the K-l-A phase as a function of the input merger param- 
eters, and summarize the results in ^ providing fitting 
formulae that can be used in semi-analytic modeling. In 
Jjni we discuss the implications of our study by comparing 
observed K-|-A fractions to those inferred by our lifetime 
calculations. We summarize and discuss the implications 
of these studies in ^ and conclude in ^ 

2. METHODS 

We employ the Gadget 2 code ()Springelll2005D to sim- 
ulate galaxy merger scenarios. In ^2.11 we describe the 
mergers studied in this work and the relevant physics in- 
cluded. Stellar light and dust extinction are modeled 
post-simulation using the three-dimensi onal polychro- 
matic radiative transfer code Sunrise (|Jonsso 3 120061: 
iJonsson. Groves, fc Coxl 120101, he reafter J06, JIO), with 
specific models discussed in ^2.21 We produce spatially- 
resolved spectral energy distributions (SEDs) along mul- 
tiple sight lines from 50x50 pixel cameras with each pixel 
spanning two physical kiloparsecs. To select K-l-A galax- 
ie s, w e mimic the data of optical spectroscopic surveys 
( ^2.3[) and generate spectral line catalogs for each simu- 
lation. For this purpose, we use the following equivalent 
widths: the Balmer lines (H/3, H7, Rd) = (486lA, 434lA, 
4IOIA) and the [O II] A = 3727A transition (EH). In 
Figure [T] we demonstrate our aperture and spectrum fit- 
ting procedure. 

2.1. Galaxy Simulations 

This study is one in a series based on numerical sim- 
ulations of gas-rich galaxy mergers. The calculations 
used here are essentially a subset of t hose employed by a 
number of othe r author s, f or exa mplelCox et al.N2006af) . 
Hopkins et all (l2005bl) . iNaravanan et al.l (l2010blla ). 

M 
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and others. The focus of this work differs in that we 
wish to generate and analyze mock spectral catalogs of 
merger remnants in order to quantify their observed post- 
starburst characteristics. Hence, most of the discussion 
in this section will be on the application of the radia- 
tive transfer code Sunrise (J06, JIO) to the extraction 



of realistic spectra. For a complete description of the 
physical models applied to the hydrodynamic methods in 
the context of this ong oing study, the reader is encour - 
age d to see [Sj jringcl. Di Matteo. fc Hernguistl (|2005bf ) 
and lCox et aLn : 200681. 

Gadget 2 (.Springel I2005D combines an accurate 
N-body gravity solver with a smoothed particle hy- 
drodynamics (SPH) code that is fully conservative 
(jSpringel fc Hernguistl 120021 ) to perform cosmological or 
galaxy simulations. Gravity dominates on galactic scales, 
but hydrodynamics is essential in many applications and 
serves as the basis for treating a great deal of chal- 
lenging baryon physics. Our simulations spawn star 
particles stochastically from the gas via a prescrip- 
tion designed to re produce the Schmidt-Kennicutt law 
(jSpringel fc Hernou ist 2003). Each galaxy is seeded with 
a central massive black hole that undergoes Eddington- 
limited Bondi-Hoyle accretion. These processes in- 
evitably inject energy back into the surrounding gas (su- 
pernovae, jets, winds, radiation, etc), so our simulations 
model this by therma lly coupling some of this energy to 
nearby gas particles fDi Matteo, Springel. & Hcrnquisti 
2005: Springel 2005; Springel, Di Matteo, & Hcrnquisl 
l2005bD . The feedback efficiency parameters used 
are identical to tho s e used in previo us studies (e.g. 
iHopkins et al.ll2005bt iCox et al]|2006al and subsequent 
works): the radiative efficiency is 0.1, and the feedback 
energy couples dynamically to the surrounding gas with 
an efficiency of « 5%. 

Our star-forming disk galaxies consist of dark mat- 
ter halos, stellar and gaseous disks, and steUar bulges. 
We place two disks on parabolic orbits with a variety 
of relative disk orientations. Generally speaking, the 
major mergers (two galaxies equal in mass) proceed as 
follows. As the two disks approach for the first time, 
strong non-axisymmetric forces cause the disks to form 
a transient bar, driying some m aterial into the center 
(| Barnes fc Hernquisti I1991L I1996D . How much material 
depends on the merger orientation as well as the proper- 
ties of th e original galaxy dH opkins et al. 2009b, d). For 
example, iMihos k Hernguistl (11994^ 119961 find that the 
presence of a stellar bulge can resist perturbations owing 
to this first passage. 

Ultimately, the torques produced during this close pas- 
sage transfer angular momentum from the halos' mutual 
orbit into their constituents, leading to a direct coales- 
cence of the two disks. During this phase, the perturbing 
forces are so great that they completely disrupt the disks, 
sending much of the remaining gas into the center of the 
combined dark matter potential. As the gas density in- 
creases in the central regions, star formation and central 
black hole growth proceed rapidly. 

In this study, we explore isolated binary mergers for 
which we vary galaxy mass, gas fraction, and bulge frac- 
tion, as well as their relative orientation and mass ratio. 

2.2. Radiative Transfer 

We generate mock optical observations of our galax- 
ies in post-processing for all relevant times during each 
merger simulation. Roughly speaking, knowledge of the 
global star formation rate and approximate wavelength- 
dependent dust attenuation is sufficient to produce a 
reasonable prediction of the integrated galaxy spectrum. 
However, large spectroscopic surveys (in which poststar- 
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burst galaxies are identified) typically measure only a 
portion of the galaxy light in a spatially-biased way 
(fibers, slits), and we wish to explore how K+A selection 
varies with aperture size and viewing angle. Therefore, 
we prefer a method that uses and preserves the infor- 
mation we get from the simulations about the locations 
of individual star particles and hence of distinct stellar 
populations, as well as the gas and dust. This requires 
knowing the complete spatially-resolved star formation 
history (SFH) of the system, which we can follow via the 
histories of the stochastically-spawned SPH star parti- 
cles. We then must assign each particle a single-age stel- 
lar population (SPS) SED for a cluster of that mass, age, 
and metallicity. Moreover, the wavelength-dependent at- 
tenuation cannot be reliably predicted from integrated 
quantities, as it depends on the geometry of stars, AGN, 
and dust, as well as dust grain composition and sizes, and 
scattering of light into the line of sight. Thus we calcu- 
late the dust attenuation along each line of sight using 
knowledge of the surviving gas particles' spatial distribu- 
tion and metallicity by assuming a dust-to-metals ratio. 

To these ends, we employ the Sunrise polychromatic 
Monte Carlo radiative transfer code (described in JOG 
and JIO). This code performs all of the necessary steps 
above and allows us to conveniently vary important stel- 
lar modeling parameters for our study. Below we discuss 
the principle components used to compile our final spec- 
tra. In addition, we describe the modeling assumptions 
we make when limited by the finite resolution of the sim- 
ulations. These choices were made to be reasonable for 
local merger remnants, and were set before making any 
of the conclusions in fJ2]or ^jSl In some cases f ji2.5p . we 
varied these assumptions to explore their effect on a small 
subset of our simulations. However, in this work we make 
no attempt to constrain or fit these modeling uncertain- 
ties, but rather treat them as a fixed set of inputs with 
which to test the overall picture. 

1. Normal stars. To model stars of age > 10^ years 
(10 Myr), we apply a high-optical-resolution (AA/A w 
20,000) STARBURST99 SED (Leithcrcr ct al. 1999) with 
metallicity closest to that of the star particle. We use 
th e two-se g ment power-law initial mass function (IMF) 
of iKroupal ()200lD . It is from these models that Balmer- 
series absorption enters our SEDs, since stellar clusters 
of a certain age will be dominated by A stars, originally 
defined by such strong absorption features. 

2. Young stars. For stars less t han 10 Myr ol d , SuN - 
RiSE uses the SED templates of iGroves et~all (|2008l ). 
These models combine one-dimensional shock models 
for evolution of the H II and photodissoci ation regions 
(PDRs) around massive stars and clusters (iDopit a et al.l 
[2005; Groves et al. 2008), with the Mappi ngs dust and 
photo i onization radiative transfer code ([Binette et al.l 
1981 ISutherland fc Dopital [1991 IGroves et al.1 120041 : 
Allen et al.ll2008l) ' ;o calculate the SED that emerges from 
the HII region and, optionally, the PDR, which is the re- 
gion where the gas surrounding young stars transitions 
from being fully ionized to fully molecular. These SEDs 
include nebular emission lines. An uncertain tunable 
model parameter of this component is the time-averaged 
fraction, /pdr, of each H II region that we assume to 
be surrounded by the birth cloud from which the clus- 
ter was born (the PDR). The dust in the PDR is a sig- 
nificant source of attenuation for young stars that are 



potentially the dominant contributor to the galaxy lu- 
minosity. Much observational and theoretical work is 
attempting to constrain this parameter, and a number 
of discussions about its impact and use for study i ng suc h 
simulations can be found in JIO. IGroves ct ah (2008^, 
iNaravanan et al.l ([20091 and I Younger et a l. (2009). 

Prior to running any simulations, we selected /pdr = 
0.3 as a realistic value for the majority of K-l-A galaxies at 
low redshifts. This choice gives results that are identical 
to within 10% (in line fiuxes) to when the fiducial value of 
0.2 from JIO is used. These values of /pdr, correspond to 
molecular cloud clearing timescales of ~l- 2 Myr, which 
is con sistent with observational estimates ([Groves et al.l 
120081 and references therein). As pointed out in JIO, the 
MappingsIII particles contribute weakly to the galaxy 
SEDs when the specific star formation rates are at or 
below those of typical local galaxies, which is what we 
expect for the majority of K-t-A galaxies found in local 
surveys. See ^2.5l for a discussion on the effect of /pdr = 
0.9 on K-l-A selection. 

Another source of noise tied to the treatment of young 
stars owes to the stochasticity of the star formation pre- 
scription in the simulation. Although it gives the ex- 
pected overall star formation rates for each gas particle, 
the version of Gadget we use for this study generates 
star particles using a finite table of random numbers for 
the Monte Carlo step. This causes the formation of star 
particles at a given time to be more correlated than it 
ought to be, and the emission line star formation indica- 
tors refiect this by being highly variable from snapshot 
to snapshot even though on average the star formation 
histories are correct. For all of our measurements, we 
smooth the spectral line catalogs over a 50 Myr box to 
limit the noisiness of the lightcurves owing to this effect. 
This will add a comparable amount (^ 50 Myr) of uncer- 
tainty to our measurements, but this is a small random 
effect. 

3. Dust attenuation. To attenuate the light originating 
from the source particles. Sunrise performs Monte Carlo 
radiative transfer on an adaptively refined Cartesian 
grid. For complete details, see J06 and JIO. We assume 
a dust-to-metals ratio of .4 (iDwek 1998) and Mil k y Way 
dust grain model s from iWeingartner fc Draind ([20011 ) 
with iDraine fc Lil ([2007[) updates. Using the energy in- 
put into the dust by the source particles. Sunrise com- 
putes self-consistently absorption, scattering, and subse- 
quent infrared re-emission from dust. The wavelength- 
dependent dust attenuation is calculated correctly from 
the assumed dust grain models, producing spectra for 
which the lines are attenuated realistically. 

The radiative transfer of dust attenuation by dense 
clumps requires subresolution considerations because the 
underlying hydrodynamic simulations do not have suffi- 
cient resolution to model th e ISM on all scales. Und er the 
ISM treatment that we use ([Springel fc Hernauistl [20031 , 
gas particles in Gadget should be thought of as contain- 
ing both a diffuse "hot" component and a clumpy "cold" 
component where the total masses of each are estimated, 
but we have no information about their physical distri- 
bution. In this work we use the PDR model described 
above to treat the significant attenuation by the CMC 
from which a star forms, but assume that the covering 
of young stars by GMCs that are not their birth clouds 
is negligible. In this approach, we account for the cold 



Poststarbursts from Mergers 



5 



phase (which can contain ^ 90% of the gas mass) only 
through the PDR model, which we use as the source spec- 
trum for newly-spawned (< 10 Myr) star particles. This 
model is sometimes referred to as "multi-phase ON" . We 
believe this choice is appropriate for the conditions of 
most local K-l-A galaxies, where large reservoirs of ob- 
scuring dust are not available. In H2.5\ we briefly discuss 
an alternative choice. 

4. Galaxy properties. We use [O II] as a star-formation 
rate indicator to select against poststarburst galaxies, 
so in order to capture the approximately correct metal- 
licity dependence of this line, we assume the gas disks 
have an initial central stellar and gas-phase metallicity 
of 1.5Zp, and gradien t —0.03 dexkpc"^, consistent with 
iZaritskv et aLTlfTOOl . We assume that the stars initially 
present in the disk were formed continuously over 14 Gyr 
and the bulge was formed in a single burst 14 Gyr ago. 

5. AGN. This implementation of Sunrise uses the 
black hole accretion rate from the Gadget simulation to 
estimate the bolometric luminosity of the central AGN, 
L — eMc^, where we use e = 0.1. The code then inserts 
an unreddened luminosity-depe ndent broadband quasar 
SED from the templates in .Hopkins et al.l ()2007c[ ). This 
source is considered to be at the location of the central 
SMBH and its spectrum is treated in the same way as a 
star particle, meaning that it is subject to dust absorp- 
tion and scattering. While the underlying AGN tem- 
plates used for this work include emission from narrow- 
line regions (NLR), the coarse interpolation used in this 
work to create the broadband AGN spectrum makes it 
impossible to measure this component to the lines (it is 
treated as continuum). Thus the emission lines in our 
spectra solely owe to star formation, and are not con- 
taminated by emission from NLRs. 

2.3. Simulated Spectra 

For the results we present here on the lifetimes of 
nearby poststarburst galaxies, we generate m ock catalogs 
desig ned to mimic Sloan Digital Sky Survey ()York et al.l 
120001 ) data. We fix the scale of our measurements by 
placing each galaxy at z — 0.1. The result of the Sun- 
rise procedure is, in general, a set of observations of 
the galaxy merger from several different camera angles. 
We choose seven viewing angles distributed uniformly 
in solid angle. One might think of each viewing angle 
for a given snapshot as a distinct measurement. Thus 

our sample consists of A'^camcras X A^snapshots X A^simulations 

measurements. The mock data consist of the spectral 
flux density on a square grid of pixels for each camera. 
For the present work, our pixels are at most 2 kpc on a 
side, which is sufficient to resolve typical spectroscopic 
apertures at z « 0.1. 

For a reasonable balance between spectral resolution 
and computation efficiency, we produce SEDs on an 
R = A/AA sa 1000 — 1200 wavelength grid for the range 
of interest, 3000-7000A. However, owing to technical lim- 
itations at the time these calculations were made, we 
have used R « 100 for the Mappings component of the 
galaxy SED ( H2.2\ component 2), oversampling it onto 
our R1200 grid. Thus the nebular emission lines are 
wider than they should be given the spectral resolution 
of the output wavelength grid, but the total flux of each 
line is conserved. This somewhat degrades the SEDs for 



which there is a substantial emission line component and 
renders impossible a two-component fitting of spectra in 
which there are comparably strong emission and absorp- 
tion components. To avoid confusion, we stick to simple 
line strength indicators and do not try to fit both emis- 
sion and absorption components. 

We conserve the flux in the emission line component, 
so we still accurately measure the evolution of the total 
flux in the line during the ingress to the poststarburst 
phase as star formation/line emission turns off and A 
stars become dominant. During the E-l-A phase and 
egress (as the A stars dominate then fade), the low- 
resolution Mappings model is sub-dominant and our 
w R1200 sampling becomes meaningful and effective for 
measuring the B aimer absorption lines as they shrink 
and disappear. This means that we produce spectra (or 
sequences of spectra) in which the effective spectral reso- 
lution varies unrealistically in both time and wavelength. 
The K-l-A selection criteria we employ will be sensitive to 
this, because ordinary Balmer absorption line strengths 
are calculated by measuring the equivalent width without 
any emission lines present, because they have been sub- 
tracted out. However, we will apply an [O II] limit and 
use typically-observed metallicity profiles, so we expect 
any Balmer emission filling to be small for systems that 
we observe as having small [O II] emission fluxes. We 
have checked this by running several simulations with- 
out the emission line component. In these cases, our 
pipeline measures absorption equivalent widths that are 
the same to within ~ 5% for all snapshots that satisfy 
the [O II] criterion. 

We measure the simulated spectrum observed by each 
camera at every time in four ways. By doing this for all 
of our simulations, we hope to get a rough idea of how 
these different observing modes can affect the selection 
of poststarburst systems formed from a variety of merger 
scenarios. 

: (Spectrum A) Our default spectrum is measured by in- 
cluding the full radiative transfer model (elements 
1-4 in the Section [2. 2p and by only using light that 
falls within a 3 arcsecond aperture. We center the 
aperture on the peak of the optical surface bright- 
ness via centroiding and weight the pixels accord- 
ing to their intersection with a 3 arcsecond diame- 
ter circle centered on this point. For z = 0.1, this 
corresponds to a diameter of approximately 6 kpc. 
See Figure [1] for a graphical description. 

: (Spectrum B) We turn off element 3 of Section [221 the 
attenuation by diffuse dust, and use the exact same 
aperture as for measurement A. The attenuation in 
the Mappings subresolution model, item 2 in §2.21 
is still included. Note that we do not re-center the 
aperture, so cases A and B measure an identical 
line of sight. 

: (Spectrum C) We measure the integrated spectrum of 
the square camera, which is 100 kpc on a side. We 
use this as a rough way to estimate the relative 
spatial concentration of the K-l-A component to the 
stellar populations. 

: (Spectrum D) We measure the integrated light as in C, 
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but we turn off element 3 of Section 
nation by diffuse dust. 



the atten- 



We point out that these methods generahze easily. Our 
pixel- weighting scheme can make spectra using arbitrary 
spectrograph aperture shapes (slits, etc), and we scale 
the aperture consistently using the angular diameter dis- 
tance to a given redshift and telescope properties. Where 
needed, we assume a flat cosmology with ilm = 0.3 and 
f^A = 0.7. In principle, it is feasible to mimic (in all but 
spectral resolution) any realistic observing mode. 

2.4. Line Catalogs 

Following a number of observational works, we use the 
rest-frame equivalent widths {Wune), and simple compar- 
isons between them, as primary indicators of poststar- 
burst activity. This equivalent width is defined (ideally) 
as 



line 



(1) 



where Ix is the spectral flux in the line at wavelength A, 
Cx is the spectral flux in the continuum at wavelength 
A, and the integral is taken separately over each spectral 
line in the galaxy's rest frame. 

Thus, our goal is to generate catalogs of this quantity 
for each simulation snapshot, camera angle (0-6, Sec- 
tion 12. 3p . and measurement mode (A-D, Section 12. 3|) . 
For the number of simulations, viewing angles, and ob- 
serving modes we wish to test in this work, we must 
measure in total a few times 10^ individual lines, so we 
require a robust, relatively quick way to automate these 
measurements. We describe our method in this section. 

Like real spectra, we do not know ahead of time the 
continuum level Cx- So for each spectrum, we must es- 
timate Cx on the wavelength range of interest. This 
computation is a source of uncertainty for any equiva- 
lent width measurement and is especially problematic for 
the [O II] and lines because they sit in a portion of 
the spectrum which varies rapidly with wavelength, espe- 
cially when Balmer-series absorption is prevalent. This 
makes any general continuum- fitting procedure difficult. 

In an attempt to minimize this issue, we follow the 
approach of many catalogs and calculate Cx using an it- 
erative 2(7 rejection filter with a window of width 200A. 
Note that these continuum fitting parameters are some- 
what arbitrary; we chose them so that the automatic 
measurement best resembles continua one might draw 
by hand near all our lines at all times during the simu- 
lation. It should be noted that this procedure somewhat 
underpredicts the true continuum level, since it is cal- 
culating an average of the flux near the lines while the 
true continuum level is unknown. Thus our absorption 
equivalent widths are slightly smaller than a more ro- 
bust calculation might find. There are many choices of 
parameters for this procedure that give reasonable but 
different results for the continuum estimation, and we 
find that this introduces uncertainty in a relatively bal- 
anced way among the lines of interest near the 4000A 
break, at the « 10% level. 

We then convolve (smooth) the continuum with a nor- 
malized lOOA box and subtract it from the spectrum and 
convolve (smooth) the result with a normalized ISA box 



to produce the input Sx for the line measurement pro- 
cedure. For each line, we approximate Equation ([l} by 
fitting to Sx a Gaussian curve with parameters Ax, cr, 
and Acenter as the line amplitude, Gaussian width, and 
center, respectively. The equivalent width for the fitted 
curve is then 

AxaV2^ 

Wat ~ -TTji V- (2) 

^ A {'^center } 

This fitting procedure has several advantages over a 
direct application of Equation ([Ij. Most importantly, 
relative uncertainties to the fit parameters are readily es- 
timated, allowing us to let the fit proceed with minimal 
preconditions. We simply throw out line measurements 
for which the line center is not sufficiently close to the 
catalog line center, the computed Wat uncertainty is very 
large, or the Gaussian fit parameters are very highly cor- 
related. Since we value speed for performing so many 
calcul ations, we apply a very simple MCMC procedure 
from .Press et all (|1986[ ). We find that a 1000-step burn- 
in followed by a 5000-step chain more than suffices to 
estimate the parameters and their relative uncertainties 
for lines that exist. 

The pipeline was tested extensively by visually com- 
paring the spectral lines to the automatically-calculated 
EW values. For lines that are not visually apparent, the 
parameter uncertainties and/or correlations are tens of 
percent, and the center of the fitted Gaussian often does 
not match the line catalog, so we reject measurements 
that have these properties. Lines that are obviously sig- 
nificant have small residuals (< 5%), and nearby/slightly 
overlapping lines can be separated effectively. This pro- 
cedure takes approximately 1 second per spectral line on 
modern processors: this is adequately fast for our pur- 
poses. Although the uncertainties we measure are not 
real in the observational sense because our models do not 
produce measurement uncertainties, their relative mag- 
nitudes are very effective at quickly estimating the ro- 
bustness of our line measurements. We will use Wat as 
the primary line strength measurement throughout this 
paper. 

In this work, for simplicity, we do not use Ha as a 
star formation rate indicator. This line enters our model 
spectra from all of the model components mentioned in 
12.21 and the resolution limitations for the emission com- 
ponents mentioned in §[2T3]are exacerbated by AGN con- 
tributions and the nearby Nitrogen lines. 

Future work on these simulations will include the Ha 
line. It is useful as an additional observational discrimi- 
nant against dust-obscured star-forming contaminants to 
the K-l-A population. Moreover, a number of different se- 
lection criteria have been applied to select poststarburst 
galaxies in different observational samples. In principle 
our models can mimic any of these, but here we do not 
consider more than simple EW cuts. This is sufficient 
to test K-l-A formation and discuss their properties, but 
not for complete comparisons to real samples. A full ex- 
periment on the merits and results of different selection 
methods is beyond the scope of this paper, but will be 
the subject of future work. 

2.5. Modeling and Numerical Considerations 

In i)2.2l we described our default model and why we 
believe this choice to be reasonable for the galaxies con- 
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Figure 1. An overview of our selection method. Optical images are created for each of seven viewing angles for all simulation snapshots 
and scaled to z=0.1. The light falling within a chosen aperture is taken to be the galaxy's spectrum in our default model. A simple k-cr 
rejection filter is applied to estimate the continuum level, and the lines of interest are then fitted with individual Gaussians. The top left 
panel here shows the logarithm of the optical flux in g-r-i colors for one of our simulated K+A galaxies. We overlay the SDSS spectrograph 
fiber scaled appropriately: this is what our default spectrum A seeks to mimic. In the top right panel we show a zoom-in of our flt to the 
H(5 line in the resulting spectrum. In the bottom panel, we show the spectrum (black), our continuum fit (dashed red), and the residuals 
(black, arbitrary offset) after fitting for He, H5, H7, H/3, H«, and [O II]. This procedure allows us to create realistically-generated spectral 
line catalogs for a large number of simulations. 



sidered here. However, we might expect poststarburst 
selection to depend sensitively on these assumptions in 
certain cases, so in this section we discuss briefly how 
our measurements are affected by different modeling as- 
sumptions and numerical resolution. These variations 
are motivated by the large uncertainties in modeling the 
dense, star-forming, obscuring ISM. Such uncertainties 
arise because many of the processes that drive the evolu- 
tion of this medium exist on scales below our resolution 
or are computationally intractable for such kinds of sim- 
ulations. Hence we are limited to using sub-grid models 
with parameter choices constrained to some degree by 
observations. 

For example, the rapidly-varying physical conditions 
of late-stage mergers may imply that the appropriate 
/pDR varies in a complicated way in both time and space. 
Emission lines are strongly suppressed from regions with 
a high PDR fraction that may exist in extremely dense 
and massive molecular environments, leading to a situa- 
tion in which the optical spectrum of a powerful starburst 
may imply instead a starburst relic. By contrast, emer- 
gent emission lines are strongest when the PDR fraction 
is low, shortening the length of time a given selection 
criterion will define a K-f A galaxy. 

Thus a value of /pdr near 1.0 may be appropriate 



for specific galaxies, such as very gas-rich interactions 
at high redshift (e.g. INaravanan et al.ll20d9l ). We believe 
it is unrealistic for this work, since we will focus on times 
of the simulations that are not extreme starbursts. How- 
ever, in order to test the sensitivity of K-l-A selection 
to high values of /pdr, a number of mergers were run 
with /pdr = 0.9. In general, the K-l-A lifetimes of these 
simulations were longer by up to a factor of two than 
those using smaller values, with a second K-l-A phase 
occurring immediately after first pericentric passage and 
prior to final coalescence (unlike the /pdr = 0.3 case). 
This can be understood because the first passage pro- 
duces a weak starburst, but star formation continues; the 
Balmer absorption strength increases in both cases, but 
with /pdr = 0.9, the emission lines from the youngest 
stars are obscured enough to be considered quiescent by 
our [O II] emission line cut. We did not find any sim- 
ulation with either choice of /pdr where we selected as 
K-l-A a snapshot during the starburst. 

Another way to deal with obscuration by dense clumps 
of gas is to turn off the PDR model and instead treat 
the entire gas mass as existing in the diffuse ISM phase 
(sometimes referred to as "multi-phase OFF"). In this 
case. Sunrise treats the entire mass of each gas parti- 
cle as having uniform density over an entire resolution 
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element. We note that this choice does not necessarily 
maximize the attenuation of young (< 10 Myr old) stars 
versus the multi-phase ON model, though it does sys- 
tematically increase the attenuation for stars with ages 
> 10 Myr. Several simulations were undertaken with 
this choice, with the effect being to attenuate most of 
the light from young stars throughout the simulation, 
including A stars. Although this choice is unlikely to be 
justified during the relatively quiescent K-l-A phase of 
galaxies, this difference represents a significant modeling 
uncertainty. 

These alternative choices may be appropriate for ex- 
treme environments. However, the duration of the K-t-A 
phase is commonly of order 100 Myr or more, signifi- 
cantly longer than the duration of the peak of the star- 
burst where the modeling of dust attenutation is most 
uncertain. Thus, while our subresolution models are im- 
perfect at following the ISM during these extreme condi- 
tions and should be treated as significant uncertainties, 
this minimal overlap cannot change general conclusions 
made about the poststarburst properties of these simu- 
lations. However, they may affect inferences made about 
specific times during the merger (e.g. the peak of the 
starburst or AGN activity), so more detailed compar- 
isons may further constrain the sub-grid models. For 
example, the K-l-A pair fraction may rely on a combina- 
tion of /pDR and progenitor bulge properties (see ^3.2\i . 
or the dynamics of K-l-A galaxies may depend on the 
feedback physics during the final merger coalescence. 

Another important component to the model is the 
assumed ISM p ressu r e sup p ort. Here, we follow 
(among others) [Cro3 ([200l . ISoringel et a ll (l2005bD. 
iRobertson et al I ()2004 l2006f ) , and [Hopkins fc QuataertI 
(|2010f ) to applv an equation of state leading to a rel- 
atively stable ISM suggestive of local disk galaxies. 
We do not vary this model here, but recent stud- 
ies have relaxed this assumption in order to simu- 
late the clumpy, turbulen t ISM observed in disk galax- 
ies at high redshift (e.g. |Elmegreen et al.' '2007^. The 
reader is encouraged to see J 3ournaud et al. (2010) and 
iTevssier. Chapon. &: BournaudI (2010) (and references 
therein) for results on how this can affect the evolution of 
the ISM in simulations of turbulent gas-rich disk galaxies 
and mergers. 

Although we are unable to fully model the clumpy 
structure of the ISM, the above model variants serve to 
bracket the plausible extremes of attenuation. The de- 
fault set described in ^2.2\ and used throughout the rest 
of this paper is observationally well-motivated, but it is 
useful to consider the effects of these different extreme 
assumptions. In summary, if we assume a maximal frac- 
tion by which young stars are attenuated by dust in their 
birth clouds (/pdr — 1), galaxies with a high star for- 
mation rate more commonly satisfy the K-l-A criterion. 
If instead we assume that all of the dust is distributed 
more diffusely throughout the ISM ( "multi-phase OFF" ), 
then there exists enough dust in our merger remnants to 
mostly obscure the A stars near the center of the rem- 
nant for a large fraction of their lifetimes. While we do 
not have a single precise appropr iate description of the 
ISM structure, such false-K-hAs (jBrown et all I2OO90 or 
dust-obscured quiescent galaxies are observationally ex- 
tremely rare at best, so we are assured that our choices 
are reasonable for the K-|-A phase of galaxies and should 



serve to encompass their real behavior. 

In addition to the modeling parameter choices, we must 
also choose at what level to resolve t he underly i ng SPH 
and radiative transfer simulations. ICox et all (|2006bf ) 
studied this for a suite of merger simulations that were 
identical in spatial resolution and equal to or worse in 
mass resolution than those considered here, but with a 
slightly different feedback model. They showed that the 
star formation histories of the Gadget merger simula- 
tions are the same to within 5% when the particle number 
is increased by factors of 2, 4, and 10. This indicates that 
our SPH particle number is sufficient to resolve the star 
formation history of the mergers. Small differences in the 
simulated spectra are expected when varying the spatial 
resolution, because the star formation law we apply de- 
pends super-linearly on the gas density, and because the 
exact distribution of stars and gas will vary with resolu- 
tion. However, the integrated quantities that will most 
directly affect the K-l-A features are typically unchanged 
at the same level as above. 

The convergence properties of the Sunrise radiative 
transfer approach were demonstrated for the integrated 
SEDs of isolated galaxies in JIO. In the present work, we 
use similar grid refinement parameters, such as Ttoi = 1, 
and a fiducial resolution setting that leads to ss 500k grid 
cells for simulations in which the gas fractions are less 
than or similar to, for example, the Sbc simulation from 
JIO. During the parts of the merger that K-|-As could 
exist in the simulations with our choice of /pdr, gas frac- 
tions are only a few percent or less and star formation 
rates are at most a few MQyr"^. In these cases, diffuse 
dust attenuation is negligible and the uncertainty in the 
output spectrum arises principally from the underlying 
stellar population methods and not the radiative transfer 
procedures. To verify similar convergence for our pipeline 
to produce the post-merger fiber spectra, we re-ran a sim- 
ulation commonly presented in this work (40% initial gas, 
IO^^Mq baryonic mass, 1-1 merger) with different reso- 
lution settings, giving 250k and 2M grid cells for the 
same simulation snapshot quoted above. For these three 
cases, in Figure [21 we plot the viewing-angle averaged 
post-merger evolution of the mean EW of the Balmer 
absorption lines H/3, H7, and US, which we will use for 
K-l-A selection. The absolute difference in this quantity 
at our times of interest for K-l-A selection is less than 
12% between the low resolution and our fiducial setting, 
and less than 2.5% between our fiducial and the high 
resolution setting. The K-|-A lifetime we calculate in all 
three cases is the same to within 0.5%. 

In addition, the K-|-A selection criteria among differ- 
ent observations are not identical, and the method we 
employ here ( !j3.3p is not a perfect match to any one 
observational scheme. Unlike many approaches, we do 
not fit stellar templates to our spectra, and owing to our 
spectral resolution limitations, we do not attempt to dis- 
entangle the absorption component of the Balmer lines 
from the emission component. Thus our EW measure- 
ments can be thought of as essentially the sum of the 
(signed) emission line flux with the absorption line fiux, 
and di ffer in this fundamenta ,l way from those of, for ex- 
ample, [Zabludoffer^l] (jl996D . when emission is present. 
When the remnant becomes absorption dominated, as 
for poststarburst galaxies, our models capture the line 
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Figure 2. Here we demonstrate the relative invariance of our 
lifetime calculations to the resolution of the radiative transfer grid. 
We show the same simulation used for row a) of Figure|4]and row b) 
of Figure[5]using three different grid refinement settings. The black 
curve is the setting used throughout the rest of this paper, and 
we find that quadrupling the number of cells for a representative 
snapshot does not dramatically alter the line strength evolution 
during our times of interest; the gray box is where our fiducial cut 
might select this simulation as a poststarburst galaxy. The K+A 
lifetimes are the same to within a few percent (when all other 
parameters are fixed). Note that our definition of equivalent width 
has absorption as negative values. 

evolution faithfully. Future enhancements to our input 
spectral models will be used to improve upon these mea- 
surements. 

Another future improvement will be to include in the 
merger models a realistic treatment of gas recycling dur- 
ing stellar evolution. In the simulations presented here, a 
fixed fraction of gas (0.1, representing ejecta from Type 
II supernova) is returned instantaneously to the ISM as 
part of the stochastic mechanism used to form stars. 
However, the complexity and short timescales of gas in- 
fall in a merger imply that the evolution of star formation 
during and immediately after coalescence will be sensi- 
tive to the amount, rate, and temperature of gas returned 
to the ISM from stars formed during the burst. More- 
over, a significant fraction of stellar mass is returned the 
ISM on long timescales via stellar winds, further motivat- 
ing the use of accurate gas-recycling methods in studying 
merger remnants such as K+A galaxies. 

3. INITIAL STUDIES 

In this section we apply the above method to a vari- 
ety of galaxy merger simulations to study their general 
properties. Our systematic study of merger parameters 
follows in 21 

In i j3.ll and Wd.'Sl we describe the basic output of our 
pipeline: the spectral line evolution of each simulation 
and the "K-|- A lifetime" . We point out in ^13. 21 that stel- 
lar bulges stabilize some merging galaxies against K+A 
creation after their first passage. In ! j3 51 and 
we employ Sunrise to study in detail the effects of view- 
ing angle, dust, and aperture size on K+A selection for 
several simulations. The influence of AGN feedback is 
discussed briefly in ^13.71 

Note that in what follows, we take the convention that 
equivalent widths in absorption are negative. 

3.1. K+A Light Curves 

The fundamental output of our procedure in Sj2] is a 
catalog of line strengths at all times during a simula- 



tion. We keep catalogs for 7 viewing angles and for our 4 
spectrum models as described in ^12.31 An illustration of 
this output is shown in Figure |3l where we compare the 
evolution of these line strengths to the global star for- 
mation history in the Gadget simulation. We plot the 
data for a single camera angle and our default Spectrum 
A. In general, the star formation histories of our merger 
simulations range between those that look like row (a) 
and those that look like row (b) . The only difference be- 
tween these two simulations is the relative orientation of 
the merging disks at the start of the simulation. Thus, as 
we vary all merger parameters, we expect to see a wide 
variety in the strength and duration of the K+A spectral 
features. 

This forms a key result that we will focus on in 21 ^ 
merger between a given pair of progenitor disks produces 
a "maximal" K+A galaxy only a fraction of the time, 
because the K+A features require a particular strength 
and shape to the star formation history. 

For mergers that do not drive a substantial starburst, 
for example row (a) of Figure [3l obvious star formation 
continues well after the nuclei merge. In the spectra we 
measure, the [O II] emission line strength does not drop 
to zero even after Gyr (and cosmological accretion may 
drive star formation at near-constant levels), and our av- 
erage of the three Balmer lines {(3, 7, 6), {H), does not 
exceed 3A in absorption. For mergers that drive a mas- 
sive burst, e.g. row (b), the gas is efficiently funneled 
to the center, where it is either consumed into stars or 
thrown out by violent feedback. This sequence of events 
leads to a rapid decline in the [O II] and Balmer emis- 
sion line strength, with a nearly simultaneous increase in 
the strength of the Balmer absorption features as the A 
stars become the luminosity-weighted dominant stellar 
component. 

For comparison, in row (c), we artificially truncate 
an isolated star-forming disk (the same disks merged 
in the first two r ows). As has been p ointed out by 
other authors fe.g. IPoggianti et al.lll999l and references 
therein), this demonstrates that a K+A galaxy as we 
have defined it is not necessarily a post-starburst sys- 
tem, but can instead be formed by rapidly shutting off 
the star formation starting from a specific star forma- 
tion rate of w lO^^^yr"^ (a few M0yr~^ for the galaxy 
we tested). Other metrics, such as colors, have been 
used to distinguish the two scenarios (ICouch fc SharplesI 
I1987I) . but we do not consider ph otometry in t h is work . 
The reader is enc o urage d to see iWuvts et al.l ()2009aD , 
iScannapieco et all ()2010[ ). and JIO for detailed discus- 
sions of common observables measured from merger sim- 
ulations using Sunrise. 

3.2. Stellar Bulges 

An unexpected property of these line catalogs is the 
strong effect produced by a stellar bulge. We naively 
expected that a simulation leading to a major starburst 
might pass through a single (or no) K+A phase. How- 
ever, our simulations of mergers between disks with no 
stellar bulge component experience two significant phases 
of enhanced Balmer absorption strength along with re- 
duced emission lines. The first occurs immediately after 
the first close passage of the two progenitors, and the 
second occurs where we expected it: after the two nuclei 
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Figure 3. A summary of different outcomes. Left: Tfie evolution of the K+A parameters as a function of time for three representative 
simulations. The dashed curve corresponds to the [O II] equivalent width measurement, and the solid curve indicates the <H> measurement, 
the average of the Balmer lines. Emission is positive in these figures. The gray box is where our fiducial cut (-5.5A) selects these simulations 
as poststarburst galaxies. Right: The global star formation rate as a function of time corresponding to the lightcurves in the left column. 
Row (a) A merger that does not induce a major starburst/shutoff scenario. Row (b) A powerful merger-induced starburst leading to 
exhaustion of the gas supply. Row (c) An isolated disk with fbulge = 0.25, fgas ~ 0.4, Mf^^ryons ~ 5 X IO^'^Mq, for which we have 
artificially shut off its star formation at t = 0. Scenarios (b) and (c) evolve through a K+A phase. 



merge and produce a s tarbur st. 

Mih os fc HernauistI Hl994bl I1996D described the effect 
of a stellar bulge stabilizing the center of a gaseous disk 
against infall induced during close passages. Since our 
default measurement mode is to observe the central few 
kpc of the brightest galaxy in the image, we are particu- 
larly sensitive to differences in the star formation history 
of the central regions of our interacting galaxies. With- 
out a bulge, a close passage of a gas disk to a major com- 
panion is enough to incite a transient bar and subsequent 
burst-quench cycle in its center, even though the galaxy 
continues to form stars at > 3Moyr~^. The bulge sup- 
presses these perturbations and preserves more gas for 
star formation prior to the final burst. The torques ex- 
perienced during the final coalescence of a major merger 



are sufficient to destroy the disk regardless of the bulge. 

This effect is shown in Figure 21 where we merge two 
identical galaxies that have a moderate stellar bulge. In 
this case, we preserve the total mass (dark matter plus 
baryons), but add a stellar bulge equal to 25% of the 
resulting baryon mass. The dynamics of the two scenar- 
ios (different only in the presence of the bulge compo- 
nent) are similar enough that the nuclear gas infall and 
starburst at the final merger lead to very similar K+A 
lightcurves after coalescence. Thus, the K-l-A history of 
a merger between two disks with significant bulges can 
be approximated by ignoring the pre-merger lightcurves 
of a simulation between two bulgeless disks. This is true 
in all cases where we have mergers both with and with- 
out bulges (about 10% of our sample), and we will make 
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use of this result for the sake of efheiency. We note that 
without bulges, the burst induced by this first passage is 
highly dependent on the choice of initial disk orientation 
(see SH). 

Although K+A galaxies have a slightly higher proba- 
bility of having a visible compan ion galaxy (~ 10% vs . 
5% for non-K+A field galaxies; lYamauchi et al.l |2008() 
than ordinary field galaxies, these bulgeless simulations 
suggest that if all merger progenitors lack a dynamically 
important stellar bulge, then the companion fraction for 
K+As should be of order 20-50%. However, a fraction 
this high is not observed, consistent with the fact that 
spiral galaxies in the nearby Universe typically contain a 
bulge. 

This companion fraction may increase at earlier times, 
where potential progenitors may have smaller bulge frac- 
tions. In addition, this prevalence of a K-|-A phase 
for our late-type mergers between first passage and co- 
alescence supports the assertion that some poststar- 
burst-|-companion systems are dynamically interacting 
and not all such systems are a superposition of an in- 
terloper with a merger remnant. 

3.3. Selection Criteria and K+A Lifetimes 

In this paper we adopt K-l-A selection criteria in the 
[O II] - (H) plane, where (H) is defined as the mean 
equivalent width of 11/3, H7, and H^. This i s similar to 
the cuts used in, for example. I Zabludoff et al.l |1996.). We 
parametrize our K-|-A lightcurves by the amount of time 
that we see a given simulation satisfying a particular set 
of selection cuts. For the K-l-A lifetimes presented in H3.4I 
and beyond, we require W[oii] < 2.0A in emission, and 
W<H> < in absorption, where x is an independent 
variable. In what follows we use —7.5 < x < —3.5, where 
negative values signify absorption. We will refer often to 
our fiducial cut, where we set x = —5.5. 

3.4. Viewing Angle 

A feature of the Sunrise method that we employ is 
the ability to vary the viewing angle of our mock ob- 
servations. Any one of our simulations may correspond 
to different observed galaxies owing to differing three- 
dimensional spatial distributions of dust, stellar popula- 
tions, and AGN relative to the viewing angle. 

The placement of our fiber depends solely on the op- 
tical flux, so changes in the dust distribution among dif- 
fering lines of sight will shift the peak of this flux around 
the image as we vary the viewing angle. Thus the stel- 
lar populations intersected by the aperture vary, and the 
extent to which we see the galaxy as characteristic of a 
post-starburst will differ. In addition, when more gas 
is present, star formation may continue at higher levels 
in various parts of the merger remnant. By looking at 
images of the merging galaxies, it is clear that in some 
cases, our aperture intersects the outer regions (for ex- 
ample, a star-forming spiral arm) of one progenitor in 
front of or behind the center of the one on which we in- 
tended to put down the fiber. In other cases, the final 
burst simply does not exhaust the gas supply sufficiently 
rapidly, so star forming clumps continue to rain into the 
center well after the peak of the burst. 

In the first column of Figure [5l we quantify this effect 
on the K-t-A lifetimes for two representative simulations. 



as well as our artificially truncated disk for comparison, 
as viewed through a typical aperture (our default Spec- 
trum A). Recall that the different spectra A, B, C, and 
D are listed in Section 12.31 and are calculated simulta- 
neously for every simulation. The uncertainty owing to 
viewing angle correlates most tightly to the gas mass 
of the disk. For our major mergers most typical of lo- 
cal spirals (initial fgas ^ 40%, or at the time of merger 
< 15%), the dispersion owing to angle is of order 10^ 
years or less. However, when the initial gas fraction is 
increased to 80%, this dispersion becomes of order 10^ 
years. In Section |4] and beyond, we will demonstrate this 
effect for a wider sample of mergers by plotting this as a 
measurement uncertainty. 

We also consider the isolated disk with its star forma- 
tion artificially cut off to zero when SFR sa 6MQyr~^: 
row (c) of Figure [5l In this case, the post-starburst 
component of the galaxy is strongly obscured when look- 
ing along the plane of the galaxy, yet easily viewable 
from higher angles. Here, dust obscuration prevents the 
younger A star population from dominating the spec- 
trum. 

In the second column of Figure [SJ we include all the 
light from the system (Spectrum C) to demonstrate that 
the magnitude of this dispersion owing to dust obscura- 
tion is strongly dependent on aperture size. This empha- 
sizes the fact that the dust in these remnants is centrally 
concentrated owing to the global gas infall experienced 
during the interaction. We describe the effect of dust in 
more detail in §3.51 and further discuss discuss aperture 
bias in several simulations in 

We conclude that when considering merger remnants 
with small amounts of residual gas, viewing angle effects 
will not introduce significant errors to the inferred stellar 
populations. However, care must be taken when making 
conclusions about systems that potentially retain a sig- 
nificant gas reservoir. 

3.5. The Effects of Dust Attenuation 

Here we briefly summarize the effect of dust attenua- 
tion on our K-)-A selection within the framework of our 
chosen subresolution models. See ^ 32. 21 and W2.h\ for a 
more thorough discussion of the underlying ISM model- 
ing and the treatment of dust in dense regions. Given 
that we track four separate spectra in each of our sim- 
ulations (i i2.3|) in which we vary both the aperture size 
and presence of dust in the diffuse ISM, it is possible to 
make approximate conclusions regarding the impact and 
distribution of dust. 

First off, note that dust influences most aspects of this 
study. When signiflcant gas and dust is present, dust 
along a particular line of sight strongly affects the visible 
K-l-A lifetime. The pre-coalescence binary K-f A phase 
( ^3.2p occurs in systems where a gas disk is still intact, so 
the presence of such a phase is viewer-dependent, and the 
viewing angle dependence of ^3.4l depends strongly on the 
gas (hence dust) content and central location. The visi- 
bility and slope of 2-D profiles in the stellar populations 
(13. 6|) will depend on the dust content and distribution. 
We argue in §3.71 that any poststarburst dependence on 
AGN feedback owes primarily to expulsion of the gas and 
dust from the central few kpc. 

To see this more explicitly, we plot in the rightmost 
column of Figure [5] the same simulations described in 
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Figure 4. Demonstration of the effect of bulges on K+A dynamics. Left: The evolution of the K+A parameters as a function of time 
for two merger simulations identical except for the presence of a stellar bulge. The gray box is where our fiducial cut might select these 
simulations as poststarburst galaxies. The simulation was chosen to be of the "strong burst" type as shown in Figure [S] and we use our 
default Spectrum A as described in the text. Right: The global star formation rate as a function of time corresponding to the lightcurves 
in the left column. Top: ftulge = 0-0. Bottom: fbulgtz = 0.25. The bulgeless merger is a K+A galaxy for more than twice as long as the 
merger with bulges, inducing a "binary K+A" phase after first passage, at t — tmerger ~ —0.7 Gyr. Despite the global star formation rate 
being 6MQyr~^ during this period, the nuclei of the galaxies experience an inflow/burst /quench cycle that happens to dominate the 
spectrograp h aperture. When bulges are present, this inflow is suppressed, yet the final starburst and K+A phase are roughly identical; 
see Mihos fc HernquistI Ill994bl . [1996) for a discussion of this phenomenon. In our simulations, this behavior of a bulge suppressing the 
flrst-passage K+A phase while leaving the post-coalescence phase intact is representative of all mergers that experience a significant burst 
owing to the first passage of the two disks. 



the previous section, except we now ignore attenuation 
by diffuse dust along each hne of sight (Spectrum B). For 
simpUcity, we focus on this single set of simulations, but 
the same effect exists (though not shown in this work) for 
all studies of SI In case (a) of Figure El the merger with 
a high starting gas fraction of 0.8, the K+A lifetimes are 
in general a few tens of percent longer for each viewing 
angle when diffuse dust is included (first column) versus 
ignored (third column). In these cases dust serves to ob- 
scure [O II]-emitting star forming regions that otherwise 
disqualify the system as a poststarburst, extending the 
K+A lifetimes. In addition, we confirm that dust is re- 
sponsible for the very large line of sight scatter in K+A 
lifetime seen in Spectrum A for this case. Note that the 
two separate groups in Spectrum B result from a partic- 
ular feature in this merger: the true instrinsic scatters of 
each group are also much smaller than the first column. 

However, in case (b) where the gas fraction begins at 
0.4, the post-burst phases do not coexist with massive 
star-forming clumps, so the K+A phase cannot be dra- 
matically extended by dust obscuration of these regions 
as in case (a). Essentially, the [O II] cut is satisfied with 
or without the presence of dust. Instead, the primary ob- 
servable effect is that the diffuse dust obscures the A star 
population after the final burst and shrinks the K+A life- 
times by up to 10% for each line of sight. Note that this 



effect may also be present in case (a), but is overwhelmed 
by the effect of star-forming clumps intersecting the aper- 
ture. These tests illustrate that the magnitude and na- 
ture of the effect of dust in the diffuse ISM on K+A 
selection depends in detail on the three-dimensional dis- 
tribution of the gas and dust with respect to the A star 
population and star-forming regions. 

3.6. Radial Profiles and Aperture Bias 



In ^13.41 we showed that K+A lifetimes are shorter 
when a larger extent of the galaxy is subtended by the 
aperture, because the fraction of the luminosity falling on 
the aperture owing to A stars, which are more centrally 
concentrated than the older stars, is reduced. Thus, gen- 
erally speaking, all of our simulations are consistent with 
K+As being the result of a starburst concentrated in the 
central few kpc. 

In the top panel of Figure [5] we demonstrate this ef- 
fect in more detail than is possible with the four types 
of spectra we will discuss in Section S) We compute the 
spectrum for a wide range of spectrograph fiber aper- 
tures, and show the spatial variation in the poststarburst 
spectrum. The first (top panel) mimics the effect either 
of using instruments with differing fiber sizes or using 
the same fiber at a different distance. The second (bot- 
tom panel) mimics the data from multiple pixels available 
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Figure 5. Demonstration of the effect of viewing angle and dust attenutation on K+A selection, emphasizing the fact that the dust in 
these remnants is centrally concentrated. We measure the length of time we view a simulation as a K+A galaxy using our default selection 
cuts. Each color represents a single viewing angle that is constant across the panels. Left: Our default Spectrum A - measurements made 
from a spectrum falling within a 3 arcsecond fiber at z = 0.1. Middle: Our Spectrum C - measurements made with the integrated light of 
the entire system viewed from different angles. Right: Our Spectrum B - measurements made ignoring the effect of attenuation by dust in 
the diffuse ISM. Each row shows one simulation; (a) gas-rich major merger leading to a strong starburst; (b) identical to (a), except with 
half the initial gas content; (c) our isolated disk with artificially truncated star formation (the same simulation as row (c) in Figure [S]!. 
Thus w e can measure the dispersion in K+A incidence with viewing angle for each of our simulations: we keep track of this scatter in 
Figures lSllll bv plotting its average magnitude as error bars along the top of each panel. 



from integral field observations. The radial dependence 
of the Balmer absorption is not flat in any of the simu- 
lations presented here. 

The curves in the top panel suggest that with a fixed 
spectrograph aperture, a K-l-A selection cut applied at 
z < 0.1 will identify a somewhat different class of ob- 
ject than when it is applied at z ~ 0.5 or higher. The 
difference in US EW between z — 0.1 and 0.5 is ~ lA 
when measuring a strong-burst merger remnant with the 
SDSS fiber. For the merger simulations (neglecting pos- 
sibly highly obscured contaminants), a constant EW cri- 
terion will typically select a more extreme population at 
higher redshifts. The effect is opposite for the artificially 
quenched disk, where an old stellar bulge dominates the 



spectrum in the very center. 

In the bottom panel of Figure [6l we demonstrate the 
ability of Sunrise to study the galaxy spectra in a spa- 
tially resolved way. We plot the spectra of seven 2 
kpc-sized pixels along a line near the nucleus of one of 
our poststarburst galaxies (the one from the solid curve, 
above). This simulation snapshot is one of our "strong 
burst" mergers and clearly shows a decline in US absorp- 
tion as we look farther away from the center. 

The gradient is steepest in the central few kpc for 
this simulated K-l-A galaxy and many of the ones pre- 
sented in this paper. Several observational studies (e.g. 
Goto et al. 2008; Pracy et al. 2009) have used integral 
field unit (IFU) spectroscopy to study the radial profiles 
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Figure 6. Top: The H5 equivalent width of single galaxies mea- 
sured through different spectrograph apertures. Here we show mea- 
surements of two major gas-rich merger remnants (solid and dotted 
curves), as well as an artificially quenched stable star- forming disk 
with f bulge = 0.25 (dashed curve). Each curve is calculated at a 
single time in the simulation by measuring the equivalent width 
of the line using a range of physical aperture sizes (listed in the 
lower X-axis), and averaging over the viewing angles that give a 
K-j-A spectrum. The upper x-axis, the galaxy's "corresponding 
redshift" , is found by fixing the angular radius of the aperture to 
1.5 arcseconds (corresponding to the SDSS fiber), and then com- 
puting the redshift for which this angle subtends the physical radius 
given by the lower x-axis. We note that the two merger remnants 
have similar trends with aperture size (signal weakens with larger 
diameter), while the truncated disk has the opposite trend at small 
radii. Bottom: We demonstrate the ability of Sunrise to make in- 
tegral field spectroscopy maps by plotting the spectrum from seven 
adjacent pixels in a line near the nucleus of the simulation studied 
with the solid black curve in the top panel. We can clearly see 
the strength of the Balmer absorption lines decreasing and older 
populations becoming dominant as we move away from the center. 



of the US absorption line, with mixed conclusions re- 
garding the presence of a gradient. While we have only 
studied several simulations in this way, our sample is con- 
sistent with a clear increase in Balmer absorption EW 
towards the center of the remnant. By studying a large 
number of such simulations, it is possible to estimate 
what the distribution of gradients is expected to be, and 
compare in detail with these sophisticated observations, 
but we leave this endeavor for a future paper. 



3.7. AGN Feedback 
Several studies, such as iWild et all 



(l2009l) and 

iBrown et al.l (|2009), have explored the connection be- 
tween poststarburst galaxies and quasar activity. While 
it is thought that AGN-induced outflows carry relatively 
little gas mass, it is not clear what the effect of this en- 
ergy deposition will be on the remaining cold gas and 
subsequent star formation. Obviously, we have chosen 
a specific model to approximate these processes, but we 



can test whether or not it predicts a difference in K-l-A 
lifetime from AGN-induced star formation quenching. 

Here we briefly demonstrate that in certain cases, the 
K-t-A feature depends weakly on the presence or ab- 
sence of AGN accretion and feedback. In the example 
we simulate here, the AGN feedback is responsible for 
dispersing the remaining gas from the central few kpc, 
reducing the optical depth to the young A star popu- 
lation, thus strengthening the Balmer absorption in the 
resulting spectrum. Dust- removal is the dominant effect: 
the burst shutdown rate and total mass are not changed 
enou gh to affect t he K+ A lines. This result is consistent 
with iWild et al.l ()2009D that at least for mergers repre- 
sentative of the low-redshift Universe, the K-l-A feature 
is not dramatically affected by AGN feedback. 

We repeat one of our commonly-used one-to-one 
merger simulations with AGN accretion turned off. This 
merger has 40% initial gas fraction, 25% bulge fraction, 
stellar mass of sa 4 x IO^^Mq, and is run on the "e" 
orbit (see below). In Figure [7] we plot the evolution 
of the Balmer equivalent width as a function of time 
through the merger, where the black solid line is our 
default setting, with AGN accretion turned on, and the 
black dashed line is with AGN accretion turned off com- 
pletely. We show a single viewing angle that is identi- 
cal in all cases and plot the global star formation rates 
for comparison. In the first panel, in which we include 
the obscuration from diffuse dust, we notice a small 
(~ lA) difference between the cases with and without 
AGN feedback, where the merger with feedback expe- 
riences stronger Balmer absorption lines. This suggests 
that AGN feedback may be contributing substantially to 
the star formation quenching. 

However, in the second panel of the figure, we see that 
the unobscured stellar populations with or without AGN 
are nearly indistinguishable after the starburst. Thus, 
the slight difference in SFR between the two cases (third 
panel) is not, by itself, dramatic enough to influence the 
K-t-A phase. Instead of directly shutting down the star 
formation, the AGN feedback serves to remove obscuring 
dust from near the young stellar population, enhancing 
the Balmer absorption features. We note that this is only 
a single viewing angle of a single merger, so further study 
is needed in order to determine whether this effect has re- 
liable observational consequences. In addition, we leave 
a systematic study of residual X-ray or LINER emission 
in simulated K-|~A galaxies to future work. 

4. SYSTEMATIC STUDY OF MERGER PARAMETERS 

In the previous section, we analyzed a number of ef- 
fects that are relevant for understanding the basic out- 
put of our radiative transfer simulations. In particular, 
we quantified the effect of stellar bulges, viewing angle, 
aperture bias, and AGN feedback. For our fixed choice of 
physical models, which we argued in ii2.5l are appropriate 
for typical local mergers, these effects are dominated by 
the physical variation in K-l-A lifetime arising from the 
variety of late-stage SFHs produced by merging a given 
pair of disks. 

In this section, we quantify this variation as a function 
of four parameters that define the initial physical setup 
of the mergers. The four parameters we vary are: (a) 
relative orientation of the disks: all else equal, this pa- 
rameter determines how effectively the merger will drive 
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Figure 7. Here we demonstrate the effect of AGN feedback on the evolution of the Balmer absorption equivalent width. The solid 
line represents our default simulation (same as in Figures |4] and |6l 40% initial gas fraction, 25% initial bulge fraction, Milky Way- mass, 
one-to-one merger) with AGN accretion turned ON. The dashed curve is an identical simulation with AGN accretion turned OFF. We find 
that the presence of AGN accretion makes almost no difference to the evolution of the post-merger spectrum when obscuration by diffuse 
dust is not considered (middle panel). However, with a realistic treatment of the surviving dust content (left panel), we see that the case 
with AGN accretion experiences a longer K-f A phase with stronger Balmer absorption than the one without this accretion. We interpret 
this as AGN feedback dispersing the surviving gas and dust, reducing the optical depth to the young stellar population in the nucleus. In 
the right panel we plot the global star formation histories of the merger. 



gas to the center of the remnant; (b) total mass of the 
progenitor galaxies; (c) gas fraction at the start of the 
merger; and (d) the mass ratio between the progenitor 
galaxies. A note about gas fractions: real galaxies will ac- 
crete matter from the IGM and replenish gas formed into 
stars, a process we do not include. Therefore, our initial 
gas fractions must be set higher than observed gas frac- 
tions for typical local galaxies. The quantity that should 
be compared to observed gas fractions is the gas fraction 
at the time of merger; this ranges between 0-30% for the 
mergers we consider, consistent with starburst galaxies 
in the local Universe. 

We point out that we have fixed the orbital parameters 
of the merger such that the progenitors follow a roughly 
parabolic orbit with a first pericentric passage at ~ 7 
kpc. Changes to this orbit will likely lead to differing 
K-t-A lifetimes and line strengths, but for simplicity we 
limit our study to this initial condition. 

We will assume for the purposes of forming the K+A 
signatures that the progenitor galaxies have a stellar 
bulge component that inhibits gas infall owing to a close 
passage, and consider only the time s after final coales- 
cence. From lYamauchi et al.l (|2008| ) we have seen that 
this correction can be as high as 5-10%, but it is within 
the uncertainties of the merger rates we will be using to 
estimate the K-l-A population. The early-type spiral ap- 
proximation is reasonable for the majority of mergers in 
the local Universe, but is likely not correct for mergers 
at high redshift or specific systems that happen to have 
low bulge fraction at present times. 

We present the full data from varying the merger pa- 
rameters separately in the four Figures [SlfTTl In each 
case, we vary a single one of our four parameters while the 
other three parameters are held fixed. When a parameter 
is held fixed, we use the following values: (a) orientation 
"e" leading to a strong burst; (b) Mbaryons = 4x lO^^M©; 
(c) /gas,o = 0.4; and (d) equal mass mergers. We sum- 
marize these trends and provide fitting functions in Sec- 
tion O For clarity, we focus only on our default spectrum 
model (Spectrum A). 

As we have seen, the intrinsic scatter owing to viewing 
angle may be large: we plot the standard deviation of 



K-l-A lifetime owing to viewing angle, averaged over all 
simulations, as a function of {H) at the top of each figure, 
and find that it can be as high as 100 Myr. This scatter 
is a complicated function of both the selection cut and 
merger scenario. For mergers that produce shorter-lived 
K-|-As (< 0.3 Gyr), the scatter is roughly a constant frac- 
tion of the K-l-A lifetime, and this can be understood 
as the effect of dust attenuating the lines by a multi- 
plicative factor. The scatter for these short-lived K+A 
galaxies tends to be larger than for the long-lived ones 
because the magnitude of the viewing-angle scatter is di- 
rectly related to the amount of gas (dust) in the system 
(Figure O, and for cases considered here, mergers that 
are efhcient at consuming gas tend to make longer-lived 
K-|-As. We point out that on average, the viewing-angle 
scatter is relatively low (less than ~0.05 Gyr in all cases) 
when considering galaxies selected by our fiducial cut of 
{H) < -5.5. 

4.1. Disk Orientation 

We assume that the relative orientation of disk galax- 
ies in the Universe on large scales is essentially a random 
variable, as is the viewing angle for which we observe 
any given galaxy: the merger rates are not fundamen- 
tally different for any subset of these tests, and so our 
simulations here characterize the scatter in K-l-A lifetime 
for a given merger event. However, the merger rate does 
depend on total mass, mass ratio, and gas fraction (pa- 
rameters b-c), so ultimately we will use the results from 
the orientation study to constrain the mean K-l-A inci- 
dence for the remaining parameters. 

As shown in Figure [8l mergers between otherwise iden- 
tical disks can produce a very wide range of K-|-A life- 
times, with a spread of over 500 Myr at our fiducial 
Balmer cut. These are the simulations in which we vary 
only the relative orientation of the progenitor spin vec- 
tors. These orientations are labeled as letters b-p at the 
top of Figure [H and the exa c t para meters can be found 
in Table 1 from iCox et aTl (|2006al ). Orientations b-h 
are commonly found in the literature and include spe- 
cial ones like the "progr ade-prog r ade" h merger, while 
orientations i-p are from Barnes] (|1992D and chosen so 
to be unbiased initial conditions. These vectors are ex- 
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Figure 8. We plot K+A lifetime measurements using Spectrum A 
(the default) as described in the text, demonstrating the dramatic 
effect that different orbital configurations have on the duration of 
the ensuing K+A phase. Here we vary the initial disk orientation of 
the galaxies, while keeping fixed the mass of the progenitor galaxies 
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mass ratio (1-1 mergers). We assume that all progenitor galaxies 
for major mergers have bulges that inhibit the binary K+A phase 
prior to final coalescence. Relaxing this assumption would lead to 
a set of simulations with an even larger spread in K+A lifetimes 
owing to the orbit-dependent nature of the burst after first passage. 
The havs along the top represent the standard deviation in K+A 
lifetime (total length = 2a) owing to viewing angle, averaged over 
this set of simulations, as shown by Figure[5]for single cases. Note: 
these bars do not represent the scatter between the curves plotted 
here. 
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Figure 9. Similar to Figure |8] except we now fix the orientation 
(e) and initial gas fraction, but vary the total mass of the merging 
galaxies. The bars along the top represent the standard deviation 
in K+A lifetime (total length = 2a) owing to viewing angle, aver- 
aged over this set of simulations, as shown by Figure [5] for single 
cases. Note: these bars do not represent the scatter between the 
curves plotted here. 
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Figure 10. Similar to Figure [8l except we now fix the orientation 
(0 for all and k for one case) and total mass, but vary the initial 
gas fractions of the progenitors. The bars along the top represent 
the standard deviation in K+A lifetime (total length = 2a) owing 
to viewing angle, averaged over this set of simulations, as shown 
by Figure [5] for single cases. Note: these bars do not represent the 
scatter between the curves plotted here. The fgas = 0.05 lifetime 
curve is zero for all values of (H) shown here. 

pected to be broadly sampled by real galaxies, and we 
see from the spread in the figure that both of these orbit 
subsets span nearly the full range of K+A lifetimes. We 
conclude that a given otherwise identical pair of gas-rich 
progenitors can merge to produce a wide variety of K-l-A 
lifetimes. 

In fact, the mean K+A lifetime inferred by this study 
is ~ 0.2 ±0.05 Gyr, roughly a factor of 5 smaller than the 
1 Gyr typical of canonical estimates! For mergers of the 
disks presented in Figure [51 our fiducial cuts find K-l-A 
lifetimes longer than 0.3 Gyr only about 50% of the time. 

In general, the long-lived K-|-As are ones in which 
a rapid, strong burst occurred, consuming much of its 
gas and producing an A-star dominated spectrum. The 
short-lived K-|-As experienced a relatively less violent 
burst, preserving more gas throughout the merger and 
creating a weaker signature in the remnant stellar popu- 
lation. The dust column will be higher in these galaxies 
and will tend to more strongly obscure the K+A signal. 
This gives a shortened observed lifetime and contributes 
to the separation between classes suggested by FigurejSK- 
This can also be seen by comparing Figure [TTk tolllb. 

4.2. Progenitor Mass 

To isolate the mass dependence, five simulations were 
run with a single disk orientation (the "e" orbit in Fig- 
ure |8]), and with a single initial gas fraction (0.4), but 
with a total mass spanning a factor of ^ 100. One 
complication is that in real galaxies, properties such as 
the mass, gas fraction, and bulge fraction are corre- 
lated. However, including all the relevant coupling ac- 
tually makes the results of such a study harder to inter- 
pret. By isolating these dependencies systematically in 
this 'one at a time' fashion, as we have done here, we 
hope to gain some physical insight into the real behavior 
of galaxies. 

In the total baryon mass range 5 x 10^ — 5 x 10^" Mq, 
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Figure 11. Similar to Figure [8] except we now fix the orienta- 
tion (e) and initial gas fraction, but vary the mass of one galaxy 
to demonstrate the trend of dec reasing K+A lifetime with in- 
creasing merger mass ratio ( i|4.4l l. Solid and dashed curves rep- 
resent different initial disk orientations that lead to a strong-burst 
merger when the masses are equal. Black - 1:1 merger, compan- 
ion fgas,o = 0.4. Blue - 3:1 merger, companion fgas,o = 0.4. 
Green - 3:1 merger, companion fgas,o = 0.8. Red - 10:1 merger, 
companion fgas,o = 0.4. In each case the more massive galaxy 
has a baryonic mass of 2 X IO^^Mq. Here we give all the merging 
galaxies bulges, so the tidal forcing owing to first passage is not 
sufficient to incite a K-l-A phase in any of them. Thus for simplic- 
ity we plot the entire evolution of the merger up to ^ 3 Gyr after 
first pericenter, allowing us to measure the response to multiple 
passages in the minor mergers. The bars along the top represent 
the standard deviation in K-l-A lifetime (total length = 2a) owing 
to viewingangle, averaged over this set of simulations, as shown 
by Figure [5] for single cases. Note: these bars do not represent the 
scatter between the curves plotted here. 



the strength and duration of the post-merger K-l-A sig- 
nature does not vary much (Figure |9|), if anything dechn- 
ing slowly with larger mass. Above roughly 10^^ M©, the 
post-merger K-f A lifetimes decline rapidly. This result 
may seem counter-intuitive because we have fixed the gas 
fraction and thus expect the burst fractions to be of the 
same order or larger. However, the K-|-A features are 
sensitive not only to the total burst fraction, but also 
the detailed shape of the burst star formation history. 
The star formation histories of the most massive merg- 
ers we simulate are in fact characterized by an extreme 
burst that may be as large a fraction of total mass as the 
low-mass cases. 

However, this extreme burst occurs primarily in the 
mode of a monotonic decrease of global star forma- 
tion, leading to an efficient, rapid merger. The final 
"starburst" in these cases is significant, but not distinct 
enough to form a long-lived K+A galaxy. By contrast, 
the less massive mergers experience a burst that better 
resembles what we might expect for making K+A galax- 
ies: an enhancement followed by a rapid decline in star 
formation. 

We note that the trend of K+A lifetime at high masses 
does not strongly affect predictions made about the pop- 
ulations observed in large surveys, which will be domi- 
nated by significantly more populous lower-mass galax- 
ies. 



4.3. Gas Fraction 

The trend as a function of gas fraction is shown in 
Figure [101 At fgas,o = 0.8, the mean K+A lifetime 
as seen by our aperture is ~ 100 Myr shorter than for 
fgas.o = 0.4. However, here we do not simulate many 
mergers with gas fractions as high as 0.8, so it is not 
clear if this trend will actually turn over for the galaxy 
population as a whole. We might expect this to happen, 
as extremely gas-rich merger remants ha ve been shown to 
reform substantial disks in some cases (|Robertson et al.l 
[200llHopkins et al.ll2009dD . and our result is then consis- 
tent with the picture that mergers with disky remnants 
experience shorter and less-robust K+A phases. As we 
showed in §3.41 a large range of K+A lifetimes will be 
seen when the gas fraction is this high, with star forma- 
tion often occurring at high levels after the merger, and 
high dust column densities obscuring various parts of the 
burst/post-burst stellar populations. 

When the gas fraction is lower, we recover an expected 
behavior. Since the gas reservoir at the time of merger 
defines an upper limit to the burst mass fraction, for low 
gas fractions at merger (< 0.1) the resulting starburst 
will not contribute enough to the spectrum to look like 
a strong K+A galaxy. 

4.4. Merger Mass Ratio 

Since equal-mass mergers are uncommon, we also con- 
sider a number of unequal-mass scenarios. We define the 
parameter fi as the ratio of the baryon mass in the more 
massive galaxy to the baryon mass in the less massive 
galaxy. In Figure [Til we plot the K+A lifetimes for a 
variety of merger mass ratios. In addition, we consider 
a second initial disk orientation in order to get a better 
idea at how this changes the results. 

The K+A lifetime falls off very quickly as we increase 
/i: on average, three-to-one mergers live as a strong K+A 
(EW < —5. 5 A) for < 100 Myr and ten-to-one mergers 
< 20 Myr (though note the high variance with viewing 
angle). As with equal mass mergers, different orienta- 
tions introduce a strong dispersion in K+A lifetimes. 

Mergers with a high mass ratio t end to make weaker 
starbursts (see jllernquist 1989; Hcrna uist fc Mihoill995l : 
iHopkins et aIlT2008d : .Cox et al.. 2QM), owing to the 
mass-dependence of the response to tidal forcing 
(|D'Onghia. E.. et al.|[2010[ ). and this is refiected in the 
shorter mean K+A lifetimes for these scenarios. In ad- 
dition, when a K+A signature is detected for a high- 
mass-ratio merger, it is often from a single viewing angle 
(out of seven), and for a single simulation snapshot at 
a single crossing between the galaxies (out of several). 
This leads to rather unpredictable behavior when the in- 
teraction is weak, likely in part because the remnants 
arc more disky than the major merger counterparts and 
hence the burst population experiences greater dust at- 
tenuation from certain viewing angles. Moreover, our 
method of placing the spectrograph fiber on the peak of 
the surface brightness leads to noisy measurements as 
the interactions proceed. The small bursts in these in- 
teractions create stellar populations with complex three- 
dimensional motions and luminosities that vary rapidly 
with time. 

5. SUMMARIES AND FITTING FORMULAE 
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We show a summary of the trends with the selective 
parameters mass, gas fraction, and mass ratio, averaged 
over viewing angle and disk orientation in Figure [T2] We 
will focus on our fiducial K+A cut, which is shown as 
the filled circles. 

One conclusion to draw from this systematic study is 
that mergers produce a wide variety of K+A behavior, 
even when the gas content and structure of the progenitor 
disks are not changed. The lifetime of this phase is not 
limited to a fixed single value for all mergers of the same 
two disks. In particular, only gas-rich major mergers 
10% gas fraction at merger) of specific initial disk ori- 
entations generate a strong K-|-A signature (> 5.5A EW 
in Balmer absorption) for a substantial fraction of an A 
star lifetime. We find that mergers reproduce canonical 
estimates of K-f A lifetimes of 1 Gyr (e.g. from SSPs) 
only when the dynamics of the merger lead to a burst 
that is in some sense "maximal" for producing the K-l-A 
spectral features. 

This means that the dispersion in K-l-A lifetime for a 
merger, given its progenitors, is large (of order 100% in 
Figure IS]). However, when comparing to a large popu- 
lation of galaxies at a snapshot in time, we care mostly 
about the estimate in the mean K-l-A lifetime for a given 
pair of progenitors, which is known better by a factor of 
3-5 because we can average over these 15 possible spin 
orientations. In Figure 1121 a bar representing the la 
scatter owing to disk orientation is shown in the upper 
left corner, and the corresponding (smaller) standard la 
error in our estimate of the mean value is shown to the 
right of it. 

In detail, our study in i j4. II finds a mean K-l-A lifetime 
that is roughly a factor of five lower than typical esti- 
mates of 0.5-1.5 Gyr: here, Tk+a ~ 0.2 ± 0.05 Gyr. In 
the next section we show that this rarity of strong post- 
starburst mergers is consistent to within a factor of a few 
with nearby large- volume spectroscopic surveys. 

For convenience, we provide example analytic approx- 
imations to our general trends that are plotted as the 
solid black curves in Figure fT2l 

As a function of baryonic mass M above roughly 
IO^Mq, we find that the mean K-|-A lifetimes for our 
fiducial selection cut (— 5.5A) follow an exponential: 



TMiMbaryon) = 0.2 CXp 



Mba 



10" Mq 



Gyr. (3) 



As a function of baryonic mass ratio /i up to at least 
^ w 10, they follow a power law: 



T^ifi) - 0.2 ^1-'-^ Gyr 



(4) 



Most of the simulations we present here have 
fgas{tmerger) ^ 0.35, SO our best-guess valucs are sig- 
nificantly uncertain beyond this regime. In the figure, 
we plot two options, one in which there is a turnover in 
the trend of K-l-A lifetimes with higher gas fractions, and 
one where it is nearly flat at higher gas fractions. As a 
function of gas fraction at the time of merger, the for- 
mer, which provides a better fit to our simulations and 
is plotted as a dotted line in the bottom right panel of 



Figure [T^ is approximated by a log- normal curve: 

Tf if gas) ^O.OM\og,oU{lo, a,) (5) 



0.045 



■ exp 



(logio(^) - ^o) 
2a2 



Gyr, (6) 



where la = —0.5 and a = 0.3. The latter option is plotted 
as a dashed black line in the figure and is approximated 

by 



TfiU 



0.55 4^ Gyr. 

1 H~ J gas 



(7) 



We derived these trends from a suite of simulations 
in which we vary only the parameter of interest in the 
progenitor galaxies. As we mentioned in §4.21 the trend 
with galaxy mass may be related to the fact that our 
more massive mergers have differing gas fractions at the 
time of merger coalescence. An added complication is 
that in real galaxies, these two properties are correlated, 
so we must be careful when combining these trends to 
ensure that we aren't mis-counting their effects. Thus, 
an appropriate combination of these trends may be one 
in which we take the gas fraction to be the driving pa- 
rameter over the mass, e.g.: 



Tifgas,f^)~0.2 Gyr- 



Tf(fgas) T^(^) 



0.2 Gyr 0.2 Gyr' 

In the next section, we use these trends to estimate the 
K-l-A fraction given merger rates from an independent 
source. 

6. IMPLICATIONS FOR THE K-l-A FRACTION 

Here we use our lifetime calculations to make a crude 
estimate of the fraction of all galaxies that are seen as 
merger-induced K-|-As. To calculate this number, we 
combine our estimated distribution of K-|-A lifetimes 
with semi-empirical merger rates. Thus we must know 
or estimate the following: (1) The K-l-A lifetime for ev- 
ery pair of possible interacting galaxies above some mass 
(flux) limit for a comparison sample, and (2) The merger 
rate for any given pair of these galaxies at the redshift 
of interest. Note that number 2 will change when the 
prevalent progenitor galaxy properties do, but number 1 
is fixed by our calculations for a given pair of galaxies. 

6.1. Merger Rates 

For item (2), the merger rates, we use the semi- 
empirical models from iHopkins et all ()2010bl i3). These 
models calculate the merger rate per galaxy per Gyr by 
applying empirical constraints to the halo occupation dis- 
tribution (HOD) of possible merger progenitors. 

We limit our analysis to merger remnants with 
M barvons > 10"^° M(T) and use the median rates quoted 
in IHopkins et al.l ()2010b[ ) for this mass range, because 
the merger rates are significantly less certain for lower 
masses and the nearby spectroscopic samples have mass 
limits of this order. 

A crude example of our calculation, undertaken in 
more detail in what follows, is to multiply a character- 
istic K-l-A lifetime by the merger rate. In the previous 
section, we found a rou gh mean lifetime of < 0.2 Gyr 
for major mergers. The IHopkins et all (|2010bl) models 
indicate that the median major (/z < 3) merger rate for 
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Figure 12. Summary of K+A lifetime trends with merger mass, gas fraction, and mass ratio. This strong dependence on progenitor 
properties drives down the ex pect ed per-merger K+A duration, leading to better ag reem ent bet ween K+A abundances and independently- 
inferred merger rates (Figure (TsJ. All points are taken from the data in Figures l9l 1101 and 1111 and the distribution of lifetimes with disk 
orientation from Figure [8[has been used to estimate the trends for the other disk pairs, assuming that the orientation is a uniform random 
variable selected from the 15 considered in this work. The filled circles correspond to our fiducial K+A selection cut {{H) < — 5.5A), and 
the curves, which are described in Section [5] represent analytic approximations to these points only. The scatter in lifetime for a specific 
merger owing to different disk orientations is > 100% (Figure [S]!, completely dominating the scatter owing to other physical effects such as 
aperture bias and viewing angle. To show this, we plot a bar representing this dominant intrinsic scatter (total length = Icr) in the upper 
left of the first panel. The smaller bar placed next to it represents the standard error in the estimate of the mean K+A lifetime averaged 
over the spin orientations obtained by dividing the larger error bar by \/l5. The size of this scatter is the same to within 30% for all three 
selection cuts. We do not plot the scatter owing to viewing angle. In the two panels on the right, the four gas fraction cases are the same. 
The only difference is that we have shifted the value of fgas to reflect the typical gas fraction at merger for each case (note: the /gas.o = 0-8 
case is sampled most sparsely in our simulations, and so both the K+A lifetime and fgas(tmerger) are more poorly constrained). 



at 



galaxies in this mass range is 0.04 galaxy ^ Gyr ^ 
z = 0.1. This implies a rough K+A incidence of 

w (0.2 Gyr)(0.04 galaxy^^Gyr^i) 
= 0.008 galaxy- ^ 



Note that our K+A lifetime per merger must still be sig- 
nificantly overestimated, because the K+A fraction in- 
ferred from large-scale spectroscopic surveys is ^ 0.001. 
This is because we have not yet taken into account the 
trends of K+A lifetime with merger parameters, e.g. 



(most importantly) the merger mass ratio. Note that 
the canonical K+A lifetime of ~ Gyr, which furthermore 
does not take into account the dispersion with orbital 
parameters, would overestimate this fraction even more 
significantly. 

6.2. K+A Fraction Estimation 

A more accurate calculation must take into account the 
trends derived in 33 because (for example), 3:1 merg- 
ers have a lower characteristic K+A lifetime than do 1:1 
mergers. The three merger parameters we consider are 
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the primary ones from the total baryonic mass Mgai 
of the merging system, the gas fraction fgas at the time 
of merger, and the baryonic mass ratio /i. 

We define the K+A hfetime to be the time the galaxy 
spend s with (H) < -5.5A, our fiducial cut, which is sim- 
ilar to the limit imposed in many observational samples. 
To estimate (1), ideally we should simulate a complete 
sample of possible mergers and apply directly our meth- 
ods to each one. In 21 '^6 performed these calculations 
for the axes of the three-dimensional merger parameter 
space that we considered, where we take the gas frac- 
tion, mass, and mass ratio as the salient properties, plus 
a range of disk orientations, for a total of w 40 simula- 
tions. To sample our entire space with a similar quality 
requires of order ~ 1000 more, so in order to reduce the 
computational expense, we limit ourselves to a rough ap- 
proximation using the information we already have. We 
expect that little more would be gained by calculating 
the complete sample. 

The general trends with each of our merger parameters 
are readily seen in Figure [T^ so we make the assumption 
that these trends extend to the unsampled sectors of the 
parameter space. For a single simulation in detail, this 
assumption is obviously incorrect. To see this, compare 
the "k" orbit from Figure |5] to the "e" orbit for 40% gas 
in Figure [TUl to the ones with 20% gas in Figure [TUl How- 
ever, for the parts of our study where we have extended 
the sample in larger numbers (e.g. Figure [TT|) . the trends 
hold to within the scatter owing to viewing angle. 

An example of this extrapolation was shown in Equa- 
tion m For the calculation here, we also include a term 
expressing the mass-dependence. However, this may not 
be appropriate in all cases, owing to the previously- 
mentioned (! j4.2p coupling between the mass-dependence 
and the gas fraction-dependence. Luckily, the steepness 
of the mass function in comparison samples renders our 
estimate of the K-l-A lifetime trend at the highest masses 
unimportant. If we assume instead a flat K-l-A lifetime 
curve with mass, thus assuming that mass and gas frac- 
tion are approximately uncorrelated, our K-l-A fraction 
estimates in what follows increase by at most ~ 10%. 

We generate a K-l-A fraction in the following manner: 

fK+A = = / P(X) T(X) i?(X) dX, (9) 

where X represents the full space of merger parameters. 
p(X) is the fraction of all merger pairs, which we esti- 
mate from observations or set by hand, represented by 
the specific choice of parameters X; this quantity has 
J p(X.) dX. = 1, and includes information about the ob- 
served merger mass function from Hopkins ct al. (2008c), 
and t he the gas fracti on distributions at z ~ 0.1 and 1.0 
from lLin et al.l ()2008[ ). 

The quantity T(X) is the K-f A lifetime derived from 
this work for the specific merger X, averaged over the 
initial disk orientation. Where we have simulations, we 
use the values from Figure [121 Where we do not have 
direct simulations, we construct T by assuming that the 
trends in Figure [12] extend trivially to the untested sce- 
narios (e.g. Equation [S]). This introduces uncertainty 
into our calculation, because we cannot be certain that 
there is not a dramatic departure from these trends. 

i?(X) is the merger rate per galaxy per gigayear from 



iHopkins et al.l ()2010b[ ). and depends on X only insofar 
as we have ignored such dependence in choosing p(X). 
We make the approximation that the universal merger 
rates at a given redshift depend only on Mgai and fi, but 
not separately on the gas fraction. Since we have chosen 
to consider a fixed mass range for which the merger rates 
are compiled, the rates we apply will depend only on /i. 

The absolute merger rates for all mergers with /i > 10 
are about twice that for "major" mergers (^ > 3) alone, 
but the K-f-A lifetimes of the major mergers are > 2 — 5 
times that for minor mergers ( §4.4^ using these defini- 
tions. Thus we might expect that the local K-l-A pop- 
ulation is composed of roughly equal numbers of major 
and minor merger remnants. 

Within the "major" merger range, our lifetime esti- 
mates are still rather different between ^J. = I and fi — 3, 
so we need an estimate of the relative size of each popula- 
tion. For this purpose, we simply choose two test values 
by hand to demonstrate the sensitivity of our models: 
iV^=3 — 3iV^=i or A''^=3 = 207V^=i. This choice of major 
merger mass ratio distribution has an effect at a factor 
of ^ 2 level, because this is the corresponding ratio in 
the K-l-A lifetime calculations. 

Since our lifetimes also depend strongly on fgas, we 
need to apply an estimate of the gas fraction distribu- 
tion in our mass range of interest. To calculate p{fgas) 
as a function of redshift, we assume gas fraction distribu- 
tions roughly consistent with observed g as-rich/gas-poor 
merger fractions from Lin et all (|2008[ ). These obser- 
vations are in good agreemen t with semi-analytic mod- 
els, e.g. iKhochfar fc IBurkertl ()2003f ). for halo masses of 
- IQI^Mq. 

For z w 0.1, where observations suggest that the to- 
tal mass density in the mass range of interest is split 
roughly equally between gas-rich disks and gas-poor el- 
lipticals, we assume that merger progenitors are 50% el- 
lipticals and 50% spirals, and strongly prefer gas frac- 
tions at merger that are < 0.1. We assume that gas-poor 
mergers do not form K-l-A galaxies (this cutoff occurs 
at roughly fgas — 0.05 for our simulations). At z w 1 
for these masses, the observed gas-free merger fraction 
is ~ 10%, and gas-rich mergers account for almost 70% 
of interactions (mixed mergers accounting for the rest). 
We assume that typical galaxies at this redshift have gas 
fractions of 0.2 by putting 50% of mergers in t o thi s bin. 

We note that the estimates from lLin et al.l ([2008) cor- 
respond to a mass range somewhat higher than the ones 
that dominate our comparison samples, so the values we 
take of the gas-rich merger fraction are probably some- 
what conservative. In addition, the choices we make of 
the gas fraction bins in which to put each type of merger 
are meant to be reasonable, but not necessarily precise. 
Changes in the exact details of the gas fraction distri- 
bution produce K-l-A fractions that vary within the un- 
certainties of other parts of our modeling. A more thor- 
ough calculation can be done using such estimates, but 
we leave this for a more precise study. 

We interpolate between our two redshift points to plot 
the K-l-A fraction as a function of redshift in the solid 
black line of Figure [T3| Such a curve should be thought 
of as the average merger-induced K-l-A fraction in all en- 
vironments at a given redshift. Note that a predicted 
merger-induced K+A fraction below the observed one 
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is not necessarily problematic, since there may be ad- 
ditional sources of K+A creation such as ram-pressure 
stripping or other quenching mechanisms. 

The gray box represents the 0.5 dex believable range 
of our best-guess model. This quant ity comes from un- 
certainty in the merger rates Hopk ins et al.l ()2010b( ) in 
conjunction with the standard error of the mean K-l-A 
Hfetimes shown in ^ We plot 0.0 < 2 < 1.2, but we do 
not attempt to account for aperture bias, and as noted 
above, our accounting for gas fraction evolution is a crude 
one (and we do not account for things like bulge mass or 
morphological sub- type) . 

6.3. Comparison to Surveys 

Above, we described our method for estimating the 
K-l-A fraction at z = 0.1 and z — 1.0. At z « 0.1, our 
procedure yields a K-l-A fraction of roughly 10"'^, and 
the K-l-A fraction we estimate at z « 1 is 0.7%, almost 
an order of magnitude higher than at z = 0.1. 

For comparison, we plot in Figure [13] the K-f A frac- 
tion with associated measurement uncertainty quoted 
from a number of spectroscopic surveys. We acknowl- 
ed ge the c omplications of this approach as pointed out 
bv lYan eFal. (2009a) and other authors, especially the 
fact that each survey uses different selection criteria, has 
different flux limits, and is affected in different ways by 
cos mic variance and enviro nmental effects. For instance, 
the IPoggianti et al.l (|1999( ) sample has a less stringent 
EW requirement than we consider here. Combining our 
models with a detailed cosmological framework, it is pos- 
sible in principle to account for such effects, but we defer 
this more sophisticated study to future work. However, 
even with this basic approach we are able to reproduce 
some expected results as a rough check of our approach. 

At z ~ 0.1, our models imply that a large fraction 
of the K-l-A population in large-sky sur veys such as the 
Las Campanas Redshift Su rvey fLCRS . iShectman et al.l 
T99I IZabludoff et al.l[l996h and SDSS (lYork et al.ll2'000r 
Gotoll2007aj ) can be accounted for by mergers. This rough 
correspondence follows directly from one phenomenon: 
that the K+A lifetimes exhibit strong dependence on 
the merger scenario. The population of mergers expe- 
riences significant variation in K-l-A incidence that de- 
pends strongly on (at least) orbit, gas fraction, and mass 
ratio. As we showed above, removing these dependen- 
cies leads to a predicted K-f-A fraction much higher than 
observed. 

Between z = 0.1 — 1.0 in the field, our example calcula- 
tion demonstrates strong redshift evolution in the K+A 
fraction. Our curve falls somewhat below the surveys 
at z > 0.2, though still within reason given the uncer- 
tainty in the models, and this might be expected if there 
is important progenitor evolution in mass, typical envi- 
ronment, bulge fraction, or if the evolution in typical gas 
fraction is more extreme than we assume. In fact, as red- 
shift increases, we expect that the progenitor galaxies in 
corresponding surveys will be somewhat higher in mass 
and later in type than lower-redshift surveys. Both of 
these effects could serve to steepen the predicted curve 
by up to a factor of 2 each toward higher redshift. An- 
other steepening effect may be aperture bias: in Figure IH 
we saw that a typical decrease in observed (H) signal be- 
tween z = 0.1 and z ~ 1 (at fixed aperture size) is « lA. 



The appropriate comparison is between the filled circles 
in Figure [12] and the squares (a stronger cut, which is re- 
quired at higher redshift to select the same galaxies) . We 
see that this could be an effect at the level of a few tens of 
percent. In addition, the other mechanisms for creating 
a K+A phase may account for some of this difference. 

7. DISCUSSION 

It is believed that the K+A criteria select galaxies 
where star formation has recently and rapidly shut down 
owing to interactions with other galaxies and/or with 
a cluster environment. The probability of becoming a 
K+A galaxy is increased if the galaxy experiences a pe- 
riod of enhanced star formation prior to the quenching, 
as opposed to the rapid truncation of a low level of con- 
tinuous star formation. Mergers and interactions are a 
natural way to create these star formation histories. We 
have studied this process through a large suite of three- 
dimensional hydrodynamic simulations coupled with the 
three-dimensional radiative transfer code Sunrise, and 
we discuss the implications of our results in this section. 

7.1. Poststarhurst Merger Models 

In combination with iWild et all (|2009| ). this work 
tested the idea that K+As arise, in many cases, via 
galaxy-galaxy mergers. Ultimately, measurements of 
K+A abundances should be reconciled with realistic 
models for their formation and event rates. The ob- 
served K+A abundance fraction among different samples 
is relatively robust insofar as the K+A selection criteria 
are uniform (though the lack of this uniformity is one of 
our uncertainties). However, realistic models for K+A 
events, specifically their duration, must be used in or- 
der to infer any information about the merger or event 
rate and hence the evolutionary context of these systems. 
We have demonstrated an order-of-magnitude calcula- 
tion that explicitly takes this into account, specifically 
calculating the K+A incidence with physical merger sce- 
nario, and we showed in 21 that the lifetime of the K+A 
phase is a strong function of both merger orbit and pro- 
genitor properties. 

In particular, the relative orientation of the spin vec- 
tors of the merging disks can be varied to create post- 
merger star formation histories with a wide variety of 
Balmer absorption strengths and lifetimes. This scat- 
ter with disk orientation represents a physical scatter in 
the mean K+A lifetime for a given pair of disks of 0.3 
dex. The K+A lifetime is maximum for equal-mass merg- 
ers and for gas fractions at the time of merger between 
^ 0.1 — 0.2, but fall of sharply at lower values. Minor 
mergers can also lead to a short, viewer-dependent K+A 
phase ( §4.4|) . While major mergers produce a more ro- 
bust, longer-lived K+A phase, the relative quantity of 
minor mergers versus major suggests that the contribu- 
tion to the merger-driven K+A population may be split 
roughly equally between major and minor. 

Thus, wc find that the average maximum K+A life- 
time for the merger of a given pair of disks is at most 
K, 0.2 ± 0.05 Gyr, a factor of at least several below the 
canonical lifetime often assumed to be ~ 1 Gyr. This is 
because not every merger event generates these canon- 
ical star formation histories. By allowing that not ev- 
ery gas-rich merger produces a K+A galaxy of order 1 
Gyr, a fixed observed K+A abundance implies a greater 
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Figure 13. An estimate of tlie K+A fraction owing to merger-induced star-formation quenching as a function of redshift, for which we 
find general agreement with similarly-selected samples. This implies that a typical merger lives as a K+A galaxy for much shorter than 
the canonical 10 years, bringing the K+A-inferred merger rate up to ones estimated by other means (,HoDkins et al.ll2010bl ). The black 
solid line is our estimate when we assume that one-to-one mergers are three times rarer than three-to-one mergers, and is approximated by 
the relation /K+A(t) oc e^t^""-^' (red dotted curve), importantly reflecting a significant increase in K-l-A abundance as z increases. The 
dashed black line is when we set this factor equal to 20; note that this choice only leads to at most a factor of 2 difference. The gray shaded 
region represents roughly the Icr confidence region of the model result, taking into account the uncertainty in the merger rates summed 
in quadrature with the scatter in the mean K+A lifetime from this work (each a factor of ~ 2 for a total ~ 3, or 0.5 dex). The points 
and error bars reflect survey data and their corresponding measurement uncertainty, with each survey using their own selection criterion. 
Error bars for the large local surveys LCRS and SDSS are lower limits to the uncertainty in the fraction computed as i/y/ Nk+a : a^nd a-U 
other error bars are as supplied by the authors. All surveys have somewhat different depth limits, corresponding crudely to lO'^" Mq of 
stellar mass, and this is the cut we use in the merger rate estimates from "Hopkins et al. (2010b). Where a survey is given two points, we 
plot both the "bright" and "faint" or total sam ple provided in the referenced paper. Our selection criteria mimic most closely the LCRS 
dZabludoff et al1ll996l) and SDSS IIGotoll2007al) survey points. A weaker cut may be more appropriate for some of the surveys for which 
an H5 cut of 3A was applied; for instance. IPoggianti et aLl II1999I ). since we include both their "k+a" and "a+k" classes in the plotted 
point. Figure [121 suggests that a factor of two correction to our model curve may be appropriate in this case (this would be insufficient). 
However, our models are much more uncertain in this regime, because a much larger class of mergers can cause such a signal, so we have 
focused primarily on ones that can satisfy our fiducial cut. 



merger rate than a calculation that assumes a lifetime 
of ^ 1 Gyr. For e xample, using a 1 Gyr duration, 
iQuintero et al.l (|2004[ ) infer a K-|-A event rate of order 
0.01 per galaxy per gigayear at z ~ 0.1, a factor of several 
below the lower limit merger rate of ~ 0.04 per galaxy 
per gigayear from lHopkins et ahl ()2010bl ) , suggesting that 
if the latter rate is correct and the typical K-l-A lifetime 
is ~ 1 Gyr, then there are too few K+A galaxies by at 
least a factor of a few. 

To reconcile the two, we have estimated the redshift 
evolution of the merger-induced K+A fraction implied 
by our s imulations in jpland F igure ll3l using the merger 
rates of iHopkins et al.n)2010bl ). For simplicity, wc con- 
sidered the strict selection cuts of lZabl udoff ct al. (19E)6|, 
e.g.), so that the observed K-l-A fraction at z = 0.1 is 
~ 10"'^. The model calculation has a large scatter of 
~ 0.5 dex, but we find broad agreement with large-area 
nearby {z < 0.2) spectroscopic samples that use simi- 
lar selection cuts, suggesting that K-f A abundances are 
consistent with a predominantly merger origin. 



7.2. Implications for Merger- Induced Evolution 

This result has significant implications for the physical 
mechanism behind the growth of bulge-dominated galax- 
ies. For example, this lower characteristic duration per 
event at fixed abundance implies that a correspondingly 
larger number of galaxies have passed through this phase. 
The mean K-f A lifetime of our simulations, weighted ap- 
proximately to the cosmological progenitor abundances 
and appropriate for this selection, is ~ 0.02 Gyr per 
event. This implies that since z = 1, '-^ 40% of galaxies 
have experienced a K+A phase, however briefly. Thus, 
many merger remnants spend a great deal of time outside 
of the traditional K-l-A spectral cuts, a consequence that 
is not too surprising, given the continuo us population of 
galaxi e s in the "p oststarburst spur" in iZabludoff et al.l 
(119961 ) : lYan et al.l (1200%^) and hints that post-burst pop- 
ulati ons may exist outside the typical sele ction regions 
(e.g. IQuintero et all 12004 lYan et al.|[2"009bD . 

Furthermore, the short-lived merger-induced K-|-As 
may be ones that experience less bulge growth than 
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the long-hved ones. Prehminary evidence suggests that 
the diskiest merger remnants in our simulations also ex- 
perience the shortest K-l-A phases, consiste nt wit h the 
picture put forth by, for example, Barnei ()2002| ) and 
[Robertso n et alJ ()20 0ff). that the merger remnant of gas- 
dominated disks will have a greater disk component. Our 
limited experiments at high gas fraction showed that 
the K-l-A phase is shorter when the gas content is large 
40%) at the time of merger, in general agreement with 
the fact that the star formation halts less suddenly, or 
not completely, in these systems. 

Hence, further study is warranted to better predict 
results for or interpret surveys of evolving galaxies at 
higher redshifts. Although gas-rich progenitors and 
mergers are more common at z > 1, we have seen that 
the K-f A incidence is a complex function of the gas con- 
tent, mass, and mass ratio of the merging pair, and the 
re-formation of gaseous disks will certainly play a larger 
role during the epoch of peak galaxy assembly at z > 2. 
Thus care must be taken to correctly interpret post- 
starburst abundance measurements in upcoming higher- 
redshift surveys. 

7.3. Resolving the Geometry and Dust 

Our combined hydrodynamics and radiative transfer 
approach allows us to study in detail the physical evolu- 
tion of galaxies and their components through the K-l-A 
phase. Roughly temporally coincident with a significant 
central starburst, SMBH accretion and feedback, and a 
rapidly-evolving ISM, merger-induced K-l-A galaxies re- 
sult from violent events that are more accurately pro- 
duced by fully three-dimensional physical models. Us- 
ing such models, we have begun to study the effects of 
dust attenuation, spatial bias, and SMBH accretion on 
poststarburst galaxies in the context of realistic mock 
observations of physically-motivated simulations. 

These studies fall into roughly two categories: one, 
understanding systematic observational effects such as 
aperture bias and viewing angle, and two, varying the 
physical setup, such as progenitor properties and mod- 
els used to evolve the structure of the ISM, e.g. forming 
stars, accreting gas onto the SMBH, and allowing for 
feedback. 

Since mergers lead to centrally-concentrated star- 
bursts, we have quantified the effect of changes in the 
spectroscopic fiber aperture between w 1 — 10 kpc on sev- 
eral simulated galaxies in ^3.6\ and shown that a change 
of order ^ lA in (H) is expected. This property may dis- 
tinguish merger-induced K-|-As from, for example, ones 
induced by ram-pressure stripping where smaller or in- 
verted gradients may be expected. 

Dust attenuation and viewing angle can introduce large 
variations in the perceived strength of K-f-A indicators, 
especially when the gas content is large. In particular, 
any diffuse dust that survives tends to obscure the central 
young stellar population and weaken the poststarburst 
signature. The scatter in the K-l-A durations induced 
by this effect is about three times smaller than the one 
owing to initial disk spin orientation 23 and the overall 
reduction in K-l-A lifetimes is an effect at the level of a 
few tens of percent (see Figure [TT|) . 

The presence of a stellar bulge suppresses the forma- 
tion of poststarburst signatures caused by relatively weak 
perturbations such as those induced during the first pas- 



sage of a merger in certain orbits (t ^3.2p . and this pic- 
ture is broadly c onsi stent w ith K-t-A pair fractions from 
lYamauchi et al.l (e.g. 120081) . This effect is dominant in a 
small aperture of 6 kpc versus the integrated spec- 
trum. Combined with knowledge of galaxy structure 
such as bulge fractions, spatially-resolved observations 
compared to simulations such as these may provide fur- 
ther constraints on the properties of merger progenitors. 

In some cases, the presence of AGN accretion and sub- 
sequent feedback strengthens the K-f-A signature of the 
merger remnant (Wild ct al. 2009, and §3.7|) . The merg- 
ers that produce strong starbursts leading to a K-l-A 
galaxy are also ones in which much gas is tunneled to 
the center of the coalescing remnant, so merger-induced 
poststarburst galaxies are likely ones that also experi- 
enced significant AGN activity. For the simulations we 
consider here, the effect of AGN feedback is limited. The 
"shutdown" we see is dominated by simple exhaustion of 
the gas supply and/or feedback related to t he sta rburst 
itself. This is consistent with iWild et al.l (|2009[ ). who 
found that reducing the feedback strength from the cen- 
tral black hole did not substantially alter the poststar- 
burst properties of such mergers. However, in at least one 
example, we find that AGN feedback produces a stronger 
K-l-A phase (^ lA EW enhancement) of longer duration 
(perhaps up to ^ 50%) than simulations without. We 
conclude that the AGN energy disperses dust and gas 
remaining in the nucleus and assists in uncovering the 
A-star population. The systematic comparison of the 
K-l-A strength distribution to a large suite of simula- 
tions may place constraints on the strength and nature 
of AGN feedback. 

Several modeling approximations were made in order 
to calculate the spectrum of simulated galaxies. These 
include a subresolution treatment for the H II and PDR 
regions surrounding young stars as well as for the state 
of the ISM within each resolution element. An impor- 
tant free parameter is /pdr, which governs the cover- 
ing fraction or clearing timescale of the molecular birth 
clouds of young stars. We chose to use /pdr — 0.3, 
corresponding to a clearing timescale of 1 — 2 Myr, as 
a reasonable approximation for local merger remnants, 
where vast quantities of molecular gas are not present. 
We have also chosen to use the multi-phase breakdown 
of the SPH gas particles to assign dust masses to either 
the hot phase or cold phase ISM. We assume that the 
attenuation of light by the cold phase ISM is captured 
entirely by the PDR model and so we do not use the cold 
phase mass to calculate additional dust attenuation. In 
^ 32. 51 we discussed several examples where we have varied 
these assumptions to see their effect, which is substantial 
in certain cases. 

7.4. Context of Present Simulations 

In work leading up to the present, similar simu- 
lations have been used to quantitatively demonstrate 
that galaxy mer gers contribut e to t he populat ions of 
z = elliptica ls (Springcl et "alll2005al: ICox et al.ll2006at 
iHopkins et aH r2006b. 2008a|7 

warm ultra-luminous in- 
frared galaxies (ULIRGs) at z « — 1 (Younger et al 
20091). dust-obscured galaxies (DOGs; i Naravanan et al 
2010bD and submillimeter galax i es ([Naravanan et al 



120091 l2010at iHavward et al] '2010'. '2011) at z > 2, and 
quasars at various redshifts (jHopkins et al..,2005a . ,2006al : 
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iHopkins fc Hernguistl [20061: iLi et all I2007D . Here, we 
have shown that K+A galaxies fit within the same frame- 
work. In particular, we have studied the idea that post- 
starburst galaxies are predominantly remnants of galaxy- 
galaxy interactions. Their properties are such that they 
will likely evolve into passive ellipticals, and these sim- 
ulations show that mergers of gas-rich disks are likely 
contributors to the K-|-A population, supporting the pic- 
ture that mergers play a significant role in driving late- 
time galaxy evolution. Thus, t he work here supports 
the gen eral picture framed by e .g. IHopkins et al.l (|2006aL 
120083 ): iSomerville et al.l (|2008| ) that these various phe- 
nomena are connected in an evolutionary sequence. 

Additional evidence for this general scenario can be 
seen in the form of a fossil record preserved in the struc- 
tural and kinematic properties of local early type galax- 
ies. For example, it has recently been shown th at the 
light p rofiles of gas-rich merger remnants ( Hopkins et al. 
2008b[) . ellipticals with central cusps (|Hopkins et al. 
2009a|), and ellipticals with cores (jHopkins et al.ll2009cD 
all show evidence of being two-component systems, with 
an inner re lic starburst left over from a gas-rich merger 
(jMihos fc Hcrnguist 1994a). These findings can ad- 
ditionally explain th e link between the yarious types 
of elliptical galaxies (|Hopkins fc Hernguistl l2010bl ). the 
contribution of merger-induced starbursts to the star 
formation history of the Universe and the infrared 
backgr ound (Hop kins et al.|[2 010a: H opkins fc Hernguistl 
I2010al: IHopkins et al.l I2006d ) . and the relationships be- 
tween superma.ssive b lack holes and their host galaxies 
(jHopkins et al.l l2007ai rbl) . The present work reinforces 
the conclusions derived in these studies by making an 
explicit connection between the immediate descendants 
of gas-rich mergers and the population of old elliptical 
galaxies. 

8. CONCLUSIONS 

We calculated rest-frame optical spectral line catalogs 
of binary galaxy merger remnants using 3-D Monte Carlo 
dust radiative transfer calculations applied to fully 3-D 
hydrodynamic simulations. Using these, we have studied 
the properties of poststarburst (K-l-A or E-l-A) galaxy 
models in an attempt to understand them in a context 
of hierarchical galaxy evolution. We summarize here our 
key findings: 

1. We tested the hypothesis that mergers of spirals 
contribute meaningfully to the population of K-|-A 
galaxies. We find that realistically-simulated ma- 
jor merger remnants commonly satisfy traditional 
K+A spectral line selection cuts O, and quan- 
tified the extent to which each of a large suite of 
merger scenarios satisfies those cuts (^J4]). 

2. The K+A lifetime is sensitive to nearly every prop- 
erty of the merging progenitors, including gas frac- 
tion, mass ratio, and especially the configuration 
of their orbit in space ( !j4.ip . This follows from the 
fact that merger dynamics exhibit a continuum of 
starburst properties, from short and violent events 
to more prolonged bursts. 

3. Typical durations of the K+A phase as commonly 
defined are < 0.1-0.3 Gyr for the merger of equal- 



mass gas-rich progenitors (SJS|), significantly shorter 
than the often-assumed value of ~ 1 Gyr. 

4. The shortness of the typical merger-induced K+A 
phase owes to the fact that not every merger or- 
bit leads to a major, succinct starburst, even when 
certain configurations of the same progenitors do 
so. Thus, the canonical model of a K+A existing 
for about 10^ years is appropriate only for the rel- 
atively few cases that the merger dynamics cause 
a sufficiently significant increase and/or sufficiently 
rapid cessation of star formation. 

5. These results reconcile measured K+A abundances 
with independently-calculated merger rates (Sj6l 
m.l^ . By assuming that all gas- rich mergers sat- 
isfy the K+A criteria for 1 Gyr, one underpredicts 
the merger/event rate by up to an order of magni- 
tude. By accurately accounting for how the K+A 
lifetime depends on physical merger scenario, as 
we have done here and for which we provide ap- 
proximate fitting formulae in ^ we find that the 
inferred merger/event rate is consistent with other 
estimates. Furthermore, our model estimates of 
the K+A abundance as a function of redshift (Fig- 
ure [13]) reflect a rapid increase from z = to z ~ 1. 

6. The subresolution treatment of dust attenuation 
remains a significant uncertainty during parts of 
the merger when much gas and dust is present 
(i i2.5|) . However, we believe our choices to be rea- 
sonable for the vast majority of K+A galaxies sim- 
ulated here, because most of the time they are rel- 
atively dust-poor. The effect of dust attenuation 
within our chosen subresolution framework is small 
and characterized throughout f|3]and 21 ^nd sum- 
marized briefly in Figure [11] and §7.31 

7. The role of AGN feedback in truncating star for- 
mation to produce the K+A features in these sim- 
ulations is minimal. The energy released by rapid 
merger-induced accretion onto the SMBH expels 
the now centrally-concentrated gas and dust from 
the merger remnant, slightly enhancing the visibil- 
ity of A stars produced by the starburst, but does 
not alter the star formation enough to be observ- 
able in the unobscured K+A signature. However, 
this work focused primarily on local analogues and 
a single feedback model, warranting further study. 

With our suite of realistic merger simulations, we 
have gained insight into how merger-induced poststar- 
burst galaxies evolve and depend on their progenitors 
and physical models. We have described and executed 
a method for deriving the true poststarburst lifetimes of 
merger remnants and reconciled K+A abundances with 
independently-estimated merger event rates. Such tech- 
niques can valuably constrain the paradigm of merger- 
induced galaxy evolution by generating observationally- 
consistent predictions for the properties of evolving 
galaxies. For example, semi-analytic modeling, based on 
approximate fitting formulae such as the ones we have 
provided here for K+A incidence, could be used to gain 
insight into related questions. Future work with sim- 
ilar techniques is warranted to disentangle the underly- 
ing causes of galaxy evolution, constrain physical models. 
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and to inform observations and theories of galaxy forma- 
tion at ever-higher redshifts. 

Specifically, the role of AGN versus star formation 
and gas consumption in the formation of ellipticals 
can be probed by varying the p hysics applie d in large 
suites of merger simulations (e.g. iWild et 311 12009') and 
comp aring them in de tail to observations of the struc- 
ture (fVang et al.ll2008f ). metal abundance fNorton et al.| 
120011 ). AGN luminosity and SEP ([B rown ct al. 2009), 
AGN/starbur st ratios (lYuan et aJ l201Qi) . and ISM 
properties (Ri ch et all 120101 ) of K-l-A galaxies, LIRGs, 
ULIRGs, and others. The complex state of the gas and 
stars in galaxies, mergers, and their remnants may de- 
pend on the precise nature of AGN and stellar feedback 
mechanisms, the details of gas recycling by stars, and the 
physics of the ISM, so that such observations and sim- 
ulations can together provide useful constraints on the 
underlying physical models. For comparisons at higher 
redshift, mergers between more gas-dominated disks with 
a turbulent ISM structure can be simulated so that we 
may interpret correctly observations of mergers and their 
remnants during the epoch of peak galaxy assembly. 
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